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Preface 


T he Bouee pour V acquisition de Series Optiques a Long Terme (BOUSSOLE) Project is a continuous ac- 
tivity, which evolves over time to deal with the multitude of activities associated with long-term oceanic 
deployments. This report is to be seen as a snapshot of the present state of the activity and achievements 1 . 
An overview of the different technical aspects is provided; it is out of the scope of this report to provide all 
details of all activities. It is directed at providing basic information that is needed for an overall understanding 
of what has been done and what is presently being done in all segments of the project. More detailed analyses 
of specific technical aspects of the project will be the subject of other, more detailed, reports. BOUSSOLE is a 
joint effort involving multiple organizations, which are combining the work of a lot of people, and is funded and 
supported by the following agencies and academic or governmental institutes: 

• Centre National de la Recherche Scientifique 2 (CNRS), 

• Institut National des Sciences de I'Univers 3 (INSU), 

• European Space Agency (ESA), 

• Centre National d’Etudes Spatiales 4 (ONES), 

• National Aeronautics and Space Administration (NASA), 

• Universite Pierre et Marie Curie 5 (UPMC), and 
• Observatoire Oceanologique de Villefranche-sur-Mer 6 . 

The resources these organizations have provided were critical in establishing BOUSSOLE as a held site for 
ocean color calibration and validation activities, but they were not the only ones who needed to cooperate. 
Indeed, a very important participant was never consulted, but luckily, decided to send a representative anyway. 
Before entering into the details of this technical report, it seems appropriate to add this silent partner to the 
introductions provided here, particularly because the task for this recently hired team member is to remain on 
site day and night to survey the mooring and provide daily status reports. After this page is turned, hopefully 
the revulsion of isolating a team member in the unforgiving environment of the open sea will be replaced with an 
appreciation for the beauty of submarine life and the largely unseen world oceanographers strive to understand. 
Incidentally, this is also recognition of the vital and hard work the technical staff associated with BOUSSOLE 
are continuously doing at sea. And, who knows, maybe you will meet this new staff member the next time 
you are at sea. This might very well reassure you — a scientist or engineer inevitably focused too much on the 
problem at hand (our common fate) — that you still possess an open mind capable of seeing what is happening 
around you and being pleasantly surprised by what is revealed as the layers of a problem are peeled away . . . 
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1 Additional details are available on the BOUSSOLE Web site http://www.obs-vlfr.fr/Boussole. 
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This is the regalec, or as it is also known, the oarfish, king 
of the herring, or ribbonfish (actually Regalecus glesne). 
The oarfish is a widely-distributed marine fish having a 
slender silvery body up to 11m (36 ft) in length, a dor- 
sal fin running the entire length of the body with red- 
tipped anterior rays rising above the head, and ventral 
fins reduced to long filaments. It swims by rhythmically 
undulating the dorsal fin (a motion resembling that of a 
snake), while keeping the body itself straight. It’s large 
size and snake- like appearance is believed to be the source 
of ancient mariner tales of sea serpents. The specimen — 
or rather, team member — shown here is about 4 m long. 
The oarfish is rarely observed in situ, but this one exhibits 
a keen interest in the status of the BOUSSOLE buoy. 


Note also the clarity of the seawater at the BOUSSOLE 
site. The details of the in-water superstructure and one 
of the support arms where the in-water instruments are 
mounted (Sect. 5.2) are clearly visible in the top-left pic- 
ture, whereas in the bottom-right picture, two of the 
divers in their small service boat and the above- water part 
of the upper superstructure are easily discerned. Oarfish, 
like this one, have been observed swimming in a verti- 
cal orientation, with their long axis perpendicular to the 
ocean surface. In this posture, the downstreaming light 
silhouettes prey, making them easier to spot. 


Photographs: David Luquet, Observatoire 
Oceanologique de Villefranche-sur-Mer 
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Abstract 

This report presents the Bouee pour V acquisition de Series Optiques a Long Terme (BOUSSOLE) project, 
the primary objectives of which are to provide a long-term time series of optical properties in support of a) 
calibration and validation activities associated with satellite ocean color missions, and b) bio-optical research in 
oceanic waters. The report starts with an introduction to the rationale for the establishing the project, including 
a definition of vicarious calibration and the specific requirements that are attached to it. The organization of 
the project is presented, and the characteristics of the measurement site, which is located in the northwestern 
Mediterranean Sea, are discussed. A qualitative overview of the collected data are provided, followed by some 
details about the new type of buoy that has been specifically designed and built for this project. Then, the 
data collection protocols, as well as the data processing techniques, are presented. The report continues with 
a quantitative summary of the collected data, as well as with a discussion of some sample results, including 
match-up analyses for the currently operational ocean color sensors, namely the Medium Resolution Imaging 
Spectrometer (MERIS), the Sea- viewing Wide Field-of-view Sensor (SeaWiFS), and the Moderate Resolution 
Imaging Spectroradiometer (MODIS). Finally, preliminary results of the vicarious radiometric calibration of 
MERIS are presented, including a tentative uncertainty budget. The results of this match-up analysis are in 
agreement with similar experiments carried out on other sites, and they allow a comparison of the performance 
of various ocean color sensors to be performed. They ultimately demonstrate the ability of the BOUSSOLE 
activity, which is combining a dedicated platform and commercial-off-the-shelf instrumentation, to provide data 
qualified to monitor the quality of ocean color products on the long term. 


1. INTRODUCTION 

A prerequisite to building a long-term (decadal) archive 
of ocean color, in response to the need for assessing the 
response of the oceanic biota to climatic changes, is to ac- 
curately calibrate the top-of- atmosphere (TOA) satellite 
observations, then to validate the surface geophysical pa- 
rameters derived from these observations, and finally to 
develop and maintain this capability over long periods of 
time. Ensuring coherence between the geophysical prod- 
ucts derived from different sensors is also an important 
aspect to consider. 

When ocean color observations from different remote 
sensors are considered as part of the task of data merg- 
ing, their intercalibration and validation might be facili- 
tated if it could be anchored on continuous long-term in 
situ stations, as recommended by the International Ocean 
Color Coordinating Group (IOCCG 1999). Deploying and 
maintaining moorings or other types of fixed stations that 
operate in a continuous way is, however, a difficult proposi- 
tion. Successful examples are provided by the Marine Op- 
tical Buoy (MOBY) program (Clark et al. 1997 and 2003), 
which has maintained an optical mooring near Hawaii since 
1995, and by the Coastal Atmosphere and Sea Time Se- 
ries (CoASTS) project maintained near Venice (Italy) since 
1995 (Berthon et al. 2002 and Zibordi et al. 2002a). 

In response to these concerns, a project was started in 
2000 with the objective of establishing a time series of opti- 
cal properties in oceanic waters to support the calibration 
and validation of ocean color satellite sensors. This activity 


has been named the Bouee pour V acquisition de Series Op- 
tiques a Long Terme] (BOUSSOLE) Project. This report 
is an attempt to provide the necessary basic information so 
that the rationale of the project and the various activities 
that are carried out are properly documented. 

1.1 General Objectives 

A new type of buoy has been specifically designed for 
the acquisition of radiometric quantities, and has been de- 
ployed in the Mediterranean Sea between the French Riv- 
iera coast and Corsica. One general objective is to perform 
match-up analyses and vicarious calibration experiments 
for ocean color satellite sensors, based on the data set that 
is being built from this new buoy. The primary objective 
of the first deployment, which took place in July 2000 and 
lasted for three months, was to validate the mooring con- 
cept. The first operational data collection occurred in May 
2002, but the entire buoy was lost because of a flaw in the 
construction of the primary flotation sphere. A redesigned 
buoy was deployed at the beginning of September 2003, 
which resulted in three months of successful data collec- 
tion. In 4 March 2004, a third data collection sequence 
was begun and lasted until July 2004, when a rotation of 
the buoy upper superstructure (i.e., the one that hosts in- 
strumentation) was carried out. Then, the data collection 
continued until June 2005, when a rotation of the entire 


f Literally translated from French as the “buoy for the acquisi- 
tion of a long-term optical series.” l ‘Boussole ” is the French 
word for “compass.” 
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buoy and mooring line was accomplished. This was actu- 
ally the beginning of a nearly uninterrupted succession of 
deployments. 

Another objective of the BOUSSOLE activity is to per- 
form vicarious calibration experiments. They should al- 
low the TOA radiance to be simulated and compared to 
the satellite measurements, in particular for the European 
Medium Resolution Imaging Spectrometer (MERIS) sen- 
sor (Rast et al. 1999). In this way, the need for a change 
in the preflight calibration coefficients for a given sen- 
sor might be evaluated, and its amount quantified. From 
this data set, match-up analyses shall be also possible for 
chlorophyll a (Chi a) concentration and water-leaving radi- 
ances, as well as algorithm evaluation (atmospheric correc- 
tion and pigment retrieval). Because of a certain common- 
ality in the band sets of the new generation of ocean color 
sensors, the data acquired with the buoy might be used for 
several of these sensors, and then contribute to the interna- 
tional effort of intercalibrating them and of cross- validating 
their products, which were amongst the basic goals of the 
Sensor Intercomparison for Marine Biological and Inter- 
disciplinary Ocean Studies (SIMBIOS) project (McClain 
1998). In addition, some protocol issues (measurements) 
are specifically linked to the use of buoys, while others, of 
general concern to marine optics measurements, may find 
specific answers when buoys are used. These aspects are 
examined, which, to the best of the project knowledge, 
have not been thoroughly investigated up to now. 

In parallel to these operational objectives (i.e., the cal- 
ibration and validation activities), the assembled data set 
will be used for more fundamental studies in marine op- 
tics and bio-optics. Among the questions and topics that 
might be addressed are the diurnal cycles of optical proper- 
ties, the response to abrupt environmental changes (storms 
and so-called red rain events), the relationships between 
chlorophyll concentration and the inherent optical proper- 
ties, the role of a fluctuating interface in establishing the 
radiative regime, the bidirectionality of the radiance held, 
the annual cycles and interannual variations of the opti- 
cal properties, the use of these properties as indicators of 
other biogeochemical properties, and the interpretation of 
the natural fluorescence signal. 

1.2 Vicarious Calibration 

The need to vicariously calibrate an ocean color satel- 
lite sensor was first demonstrated for the Coastal Zone 
Color Scanner (CZCS), which was a proof-of-concept mis- 
sion with an operational capability spanning 1979-1986. 
On one hand, this instrument was not equipped with the 
necessary onboard devices for the monitoring of the long- 
term degradation of the instrument response sensitivity 
(internal lamps suffered a rapid degradation), and, on the 
other hand, the mission did not include an extensive vi- 
carious calibration program. In the end, the calibration of 
the CZCS was never based on a robust scientific capability 
and has never been sufficiently confirmed. 


The converse situation is illustrated by the Sea- viewing 
Wide Field-of-view Sensor (SeaWiFS) mission (Hooker et 
al. 1992). The instrument is equipped with a solar dif- 
fuser and has the capability of viewing the full moon, both 
of which are used to track the temporal stability of the 
instrument; and, in addition, an extensive vicarious cali- 
bration and algorithm validation program was set up prior 
to launch (McClain et al. 1992). The vicarious calibration 
portion of the latter is based on the deployment of a perma- 
nent optical buoy, the MOBY program (Clark et al. 1997), 
which is moored in Case-1 waters near the island of Lanai 
in Hawaii. In parallel to this central and key element, ex- 
tensive campaigns are conducted around the world ocean, 
in order to collect radiometric measurements and ancillary 
data in a variety of environments (e.g., Robins et al. 1996, 
Aiken et al. 1998, and Barlow et al. 2003). These data 
are used to permanently evaluate the quality of the level -2 
products and improve, whenever possible, the applicable 
bio-optical algorithms. 

Without ground-truth data — more properly sea-truth 
data — it is impossible to maintain the calibration of a satel- 
lite sensor at the desired level over the full course of the 
mission, which is generally designed to be on the order 
of about five years (although many satellites operate for 
longer periods of time as evidenced by SeaWiFS, which 
was launched in 1997). This is not caused by any weak- 
ness in the sensor nor in the onboard calibration devices; it 
is simply because of unavoidable physical considerations: 
the goal of modern ocean color sensors is to provide the 
water-leaving radiance in the blue with a 5% accuracy over 
oligotrophic, chlorophyll-depleted, waters (Gordon 1997), 
which can be expressed as well as an uncertainty of 0.002 
in terms of reflectance (Antoine and Morel 1999). Because 
this marine signal only represents about 10% (at most) of 
the TOA radiance (i.e., the radiance directly measured by 
the spaceborne sensor), achieving this goal requires the in- 
struments involved are calibrated to better than 0.5%, or 
approximately 1%. This is extremely challenging consid- 
ering the present technology, and probably will remain an 
elusive accomplishment. 

A successful strategy for ensuring high-quality obser- 
vations from an ocean color satellite mission is based on a 
rather pragmatic approach and consists of a) making the 
best possible effort when initially characterizing and cali- 
brating the spaceborne instrument on the ground, b) using 
onboard calibration devices or maneuvers to track sensi- 
tivity changes over time while on orbit, and c) adjusting 
individual channels to force agreement with the sea-truth 
data — the basis of the vicarious calibration procedure — 
which produces a final adjustment for the whole (instru- 
ment plus algorithms) system. 

1.3 Moorings as Long-Term Data Sources 

Maintaining a permanent optical mooring is a costly, 
but pertinent, solution to the problem of how best to col- 
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lect a significant number of the needed sea-truth obser- 
vations for the vicarious calibration process. This is eas- 
ily demonstrated by comparing this option to the usual 
oceanographic practice of deploying research vessels to var- 
ious parts of the world ocean. 

Past experience has demonstrated that a one-month 
cruise is able to produce a maximum of a dozen of match- 
up points when conditions are extremely favorable. This is 
simply because of the number of conditions that must be 
simultaneously fulfilled for a measurement to be usable for 
the calibration and validation process (clear sky, low wind 
speed, acquisition within a maximum of 1 h of satellite over 
pass, nominal operation of the instruments, etc.). 

A permanent mooring is well adapted to maintaining 
a consistent time series of in situ measurements over a 
long period of time, because the needed personnel, equip- 
ment, and methodologies remain mostly fixed once they 
are established. Although the latter is a demonstrable ad- 
vantage, ensuring the same level of consistency between 
the people, equipment, and protocols from different cruises 
inevitably adds some extra uncertainties in the overall en- 
terprise of data collection and processing. A permanent 
station is also well suited for developing and testing new in- 
strumentation as well as new algorithms, and therefore, to 
permanently improve the quality and the variety of prod- 
ucts that can be derived from the ocean color observations 
at the TOA level. It is also a unique opportunity to es- 
tablish the cross calibration between different sensors by 
anchoring them to the same in situ time series. 

A criticism that is often made of the mooring option 
is the uniqueness of the measurement site. This is actu- 
ally not an argument, because the vicarious calibration 
process is a physical process, which does not require that 
data from a variety of different environments are collected. 
It requests, however, that the maximum of information 
(quantitatively as well as qualitatively) are collected at 
the selected site so that the reconstruction of the TOA 
signal through radiative transfer calculations is performed 
with the best possible accuracy. This does not preclude, 
however, that data are collected in other environments for 
validation of level-2 products and for verification of the 
calibration obtained at the initial mooring site. 

The scheme briefly exposed above (i.e., a permanent 
calibration site and a parallel program of validation from 
ships) was adopted by the SeaWiFS Project. It has been 
a success, and this sensor is probably the best-calibrated 
ocean color sensor the community has ever had, at least 
until next-generation sensors such as MERIS and MODIS 
are operational over a similar time span and demonstrate 
their full potential. 

1.4 Specific Protocol Requirements 

It has just been said that moorings are adapted to pro- 
vide large numbers of in situ data points, however, what 
would be the advantage in case these data are not of the 


desired quality ? The question arises, because the measure- 
ment conditions on a mooring may degrade the final qual- 
ity of the data. Buoys are often moving a lot, because the 
sea is rarely totally flat, so the instruments are (at some 
level) in an unstable situation. They are often made of 
a large floating body with the equipment installed under- 
neath, so the light sensors can be significantly shaded. The 
prolonged immersion of the instruments is favorable to the 
development of biofouling. Finally, the number of calibra- 
tions that can be performed, which is particularly impor- 
tant to track the temporal stability of the radiometers, is 
dependent on the frequency of maintenance visits (during 
which instruments can be exchanged). The servicing in- 
terval might be insufficient in some cases, either because it 
is inherent to the project organization or because of severe 
weather. Consequently, the qualification of a given moor- 
ing program with respect to the ocean color calibration and 
validation requirements should consider the above points 
and evaluate whether or not they have been accounted for 
through a protocol. The data collected by the moored in- 
struments should be, for instance, compared to the same 
data derived from ship-deployed instrumentation. 

2. STRATEGY 

The BOUSSOLE Project is composed of three basic 
and complementary elements: a) a monthly cruise pro- 
gram, b) a permanent optics mooring, and c) a coastal 
Aerosol Robotic Network (AERONET) station. Each of 
these three segments is designed to provide specific mea- 
surements of various parameters at different and comple- 
mentary spatial and temporal scales. When combined, 
they provide a comprehensive time series of near-surface 
(0-200 m) oceanic and atmospheric inherent optical prop- 
erties (IOPs) and apparent optical properties (AOPs), re- 
quired to accomplish several objectives: 

• Performing the vicarious radiometric calibration of 
satellite ocean color sensors (i.e., the simulation of 
the TOA radiances recorded by the sensor in var- 
ious spectral bands, including the visible and near 
infrared (NIR) domains); 

• Performing the validation of the level-2 geophysi- 
cal products that are produced from the observa- 
tions of these satellite ocean color sensors, includ- 
ing the normalized water-leaving radiances, the pig- 
ment concentrations, and the aerosol optical thick- 
ness (AOT) and types; 

• Performing both of these operations on a long-term 
basis (i.e., at least for the duration of the MERIS 
mission, for which the BOUSSOLE project has been 
primarily set up) ; and 

• Making progress in several domains of ocean optics 
and bio-optics, as already stated earlier (Sect. 1). 

A basic description of the three elements and the opera- 
tions carried out there is given below. 
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2.1 Monthly Maintenance Cruises 

From July 2001 to December 2003, the BOUSSOLE 
optics cruises were planned with a duration of three con- 
secutive days distributed on a near-monthly basis within 
the budget of 30-31 sea days per year aboard the research 
vessel (R/V) Tethys II. In these first two and a half years, 
some cruises were unsuccessful in data collection because of 
bad weather preventing ship departure or scientific work 
for the entire three days. For this reason, for 2004 the 
cruise duration during the seasons with characteristically 
unsettled conditions, October to April, were increased to 
up to five days with a subsequent reduction in the over- 
all number of cruises per year. The same logic is followed 
since 2005, yet without reduction of the number of cruises 
(i.e., 12 cruises per year). 

The so-called optics day (defined as the period dur- 
ing which the Sun is at an angle greater than 20° above 
the horizon) commences and ends with a 400 m conduc- 
tivity, temperature, and depth (CTD) profile. The CTD 
sensors are mounted on a rosette equipped with 11 Niskin 
bottles. Spectral absorption and attenuation, a(z, A) and 
c(z, A), respectively (z indicates depth and A denotes wave- 
length), is measured at nine wavelengths using a West- 
ern Environmental Technology Laboratories (WETLabs), 
Inc. (Philomath, Oregon) 25 cm pathlength AC-9 with en- 
hanced data handling (AC-9+) mounted in the twelfth bot- 
tle position. Additional sensors are the CTD and O 2 sen- 
sors, a colored dissolved organic matter (CDOM) fiuorom- 
eter, a chlorophyll fluorometer, and a backscattering (bf,) 
sensor (Sect. 6.2). 

For at least one CTD cast each day, the bottles are 
fired at selected depths during the ascent of the CTD over 
the 0-400 nr depth range. The collection depths are de- 
cided during the descent of the CTD after observation of 
the fluorescence profile, thus optimizing the representation 
of any features in the water column. For other casts, trip- 
licate samples are collected from 5 and 10 nr. In the sum- 
mer, when the optics day is longer or when there are ideal 
satellite nratch-up conditions, an additional CTD profile is 
performed around solar noon or in time with a SeaWiFS 
or MODIS overpass. 

Seawater samples (2.8 L) are filtered through What- 
man GF /F Liters using a low pressure vacuum and stored 
in liquid nitrogen. Back in the laboratory, the Liters are 
later analyzed by high performance liquid chromatography 
(HPLC) and spectrophotometry for pigments and partic- 
ulate absorption, respectively (Sect. 7.1 and 7.2). 

Multispectral upward and downward irradiance pro- 
Lles, E u (z, A) and Ed(z, A), respectively, are performed 
with the objectives of providing synchronous in situ Sea- 
WiFS and MERIS calibration and validation proLles, char- 
acterizing the light held throughout the day at the BOUS- 
SOLE site, and providing a support data set for observa- 
tions from the BOUSSOLE buoy. During these proLling 
sessions, multiple proLles are performed with a SeaWiFS 


ProLling Multichannel Radiometer (SPMR), if possible, to 
improve the quality of the data by replicates. The Lrst 
SPMR session of the day begins after the Lrst CTD pro- 
Lle, and continues ideally until the end of the optics day, 
before the Lnal CTD proLle. A gimbaled 47T photosynthet- 
ically available radiation (PAR) sensor positioned on the 
foredeck and operated from the CTD computer serves as 
a light held stability indicator during SPMR proLling. 

For the satellite passes, whenever possible, SPMR pro- 
Lling is performed within 1 h of satellite overhead passes 
of SeaWiFS and MERIS and around solar noon. Optimal 
conditions for these measurements are low humidity, blue 
skies and Lat, calm sea surface. If the sky is clear and sea 
conditions are calm, a second-generation Satellite Valida- 
tion for Marine Biology and Aerosol Determination (SIM- 
BADA) measurements are performed consecutively where 
possible with SPMR proLles (only carried out from July 
2001 to July 2003). If sea conditions are poor but the sun is 
clear of clouds, SIMBADA sun photometer measurements 
can be made at intervals throughout the day to measure 
atmospheric optical thickness. 

When sea surface conditions are calm, a low-volume 
(small shadow) surface Loat is used to suspend the SPMR 
in a vertical position with the E u sensor submerged ap- 
proximately 20 cm below the surface. The SPMR is held 
in this orientation for a duration of at least 3 min and is 
released when the light held is expected to be stable. This 
decision is a judgement call based on sky observation and 
monitoring of the PAR sensor. To start the descent, an 
electronic trigger mechanism is Ltted to the surface Loat, 
which can be controlled from inside the lab. Multiple de- 
scents ideally will be started in this way and the data will 
be used to assess near-surface upwelled radiance (L u ) ex- 
trapolation model calculations. 

For each cruise, at the end of the on-site optics mea- 
surements, there is a CTD and IOP transect between the 
BOUSSOLE site and the port of Nice (France) consisting 
of six Lxed locations. The aim is to have a representative 
proLle of the water column on either side of the bound- 
ary to the Ligurian Current. The approximate time of 
day that this transect is performed is kept similar for each 
cruise, whenever possible, to minimize the inLuence of di- 
urnal variability. 

A study of the spatial variability at the BOUSSOLE 
site was performed once on each cruise until the middle 
of 2005. A Lxed square mile quadrilateral grid based on 
Global Positioning System (GPS) data is covered by the 
ship at a speed of 7-8 kts. During transit around this grid, 
water is continuously being pumped from an inlet beneath 
the ship’s hull, and directed to the thermosalinograph and 
a Luorometer. Samples are collected at three equidistant 
points along this grid, for subsequent HPLC analyses, in 
order to convert dimensionless Luorescence measurements 
into chlorophyll concentrations. 

On other uninterrupted transits between Nice and the 
BOUSSOLE site, particularly when there is a high sun an- 
gle and clear skies, CIMEL Electronique (Paris, France) 
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Fig. 1 . The coastal area of the northwestern part of the Mediterranean Sea showing a) the southern coast 
of France and the island of Corsica, plus the generalized work area in the Ligurian Sea (black star) for the 
BOUSSOLE activity; and b) a magnification showing the position of the BOUSSOLE mooring plus the 
Dynamique des Flux Atmospheriques en Mediterranee (DYFAMED) site and the Meteo buoy maintained by 
the French weather forecasting agency (Meteo France). The positions of the six stations, which are sampled 
once a month during transits from Nice to the BOUSSOLE site, are also displayed. 


model 317 (CE-317) sun photometer measurements are 
taken at 30 min intervals (approximately 5 nmi) to charac- 
terize the variability in atmospheric optical thickness be- 
tween the Cap Ferrat sun photometer site and the BOUS- 
SOLE site. 

2.2 The Mooring Site 

The mooring site is equipped with a new type of buoy, 
specifically designed for collecting optical data (Sect. 5) 
with two specific objectives: a) obtain a near-continuous 
sampling at the surface (depths less than 10 m), in order 
to efficiently support the validation of the geophysical pa- 
rameters derived from ocean color remote-sensing obser- 
vations collected by a variety of satellite sensors; and b) 
create a high-resolution data set to support fundamental 
work in marine optics, e.g., bidirectionality of the radiance 
held emerging from the ocean, short-term changes of op- 
tical properties, response of optical properties after strong 
environmental forcings (red rains and storms), effect of 
wind-induced bubbles, near-surface behavior of radiomet- 
ric quantities, etc. 

The satellite sensors that have been or are being sup- 
ported include MERIS from the European Space Agency 
(ESA), SeaWiFS and MODIS-A (on the Aqua spacecraft) 
from the National Aeronautics and Space Administration 
(NASA), plus the third Polarization and Directionality of 
the Earth Reflectance (POLDER-3) sensor aboard the Po- 
larisation et Anisotropie des Reflectances au sommet de 
1' Atmosphere, couplees avec un Satellite d’ Observation em- 


portant un Lidar (PARASOL) satellite, from the French 
Centre National d Etudes Spatiales (ONES). 

2.3 The Coastal Site 

Since 3 July 2002, a coastal site has been equipped with 
an automatic sun photometer station, introduced within 
the AERONETf. This site is the Cap Ferrat, in front of 
the Laboratoire d'Oceanographie de Villefranche (LOV), 
by 7°19'E,43°41/N. This equipment provides a continu- 
ous record of the sky radiances (principal plane and al- 
mucantar) and of the attenuation of the direct solar beam, 
from which aerosol types and aerosol optical thickness will 
be retrieved. The annual calibration is managed by the 
AERONET group, as well as the inversion of the sun pho- 
tometer measurements in order to get the aerosol optical 
thickness. This site is established in order to provide the 
last element of data that is needed in the process of vicari- 
ously calibrating ocean color satellite observations, namely 
the aerosol types and optical thickness. 

3. MEASUREMENT SITES 

The following three subsections describe the main char- 
acteristics of a) the offshore BOUSSOLE site (32 nautical 
miles from Nice), b) the sampling done on the transect 
from the BOUSSOLE site to the Nice harbor, and c) the 
AERONET station that is installed within the military 
premises of Cap Ferrat. 

f Further details about AERONET sites are available at the 
following Web address: http : // aeronet . gsf c . nasa . gov. 
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3.1 The BOUSSOLE Site 

The site where the mooring is deployed and where the 
monthly cruises are carried out is located in the Ligurian 
Sea, one of the sub-basins of the Western Mediterranean 
Sea (Fig. 1). Water depth varies between 2,350-2,500 nr in 
this area, and it is 2,440 nr at the mooring point, which is 
located at approximately 7°54'E,43°22 / N. 

3.1.1 Bottom Depth 

Accurate bathymetry data was a critical requirement 
for properly deploying the taut mooring used to anchor 
the buoy to the ocean floor (Sect. 5). The water depth 
at the proposed site was determined by performing sev- 
eral deep CTD casts, from which the measurement of the 
pressure (translated into depth using the simultaneous ob- 
servations of temperature and salinity), was added to the 
measurement of an altimeter mounted on the CTD rosette 
(the CTD was stopped about 10 m before the bottom was 
reached). The water-column thickness was, therefore, de- 
termined with a less than 1 nr accuracy. These measure- 
ments were used to calibrate a two-dimensional mapping 
of the sea bottom that was obtained from the ship echo 
sounder. This mapping revealed a flat bottom character- 
ized by a weak northwest-to-southeast slope (Fig. 2). 
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Longitude 

Fig. 2. The bathymetry, in meters, of the area 
around the BOUSSOLE site (at the crossing of the 
horizontal and vertical lines). The open circles indi- 
cate the location of the individual deep CTD casts 
that allowed the echo-sounder survey (grey dots) to 
be calibrated in terms of water depth (gray-scale 
bar to the right of the plot). 

The tidal amplitude, the dynamic height changes, and 
the water level variations caused by atmospheric pressure 
changes are all less than about 50 cm for the deployment 



site, which is not significant for the type of mooring design 
envisioned here. They might become problematic in those 
cases wherein a similar mooring is deployed in another wa- 
ter mass where such changes are larger. 

3.1.2 Physical Conditions 

The BOUSSOLE site was selected in particular because 
the prevailing ocean currents are usually extremely weak. 
This peculiarity is due to the selected position being rather 
close to the center of the cyclonic circulation that char- 
acterizes the Ligurian Sea. The northern branch of this 
circulation is the Ligurian Current, which manifests as a 
jet flowing close to the shore in the northeast-to-southwest 
direction which, in turn, establishes a front whose posi- 
tion varies seasonally (closer to shore in the winter than in 
the summer). The southern branch of the circulation is a 
southwest-to-northeast current flowing north of the island 
of Corsica; the eastern part of the circulation is simply 
imposed by the geometry of the basin. 

The dominant winds are from the west to southwest 
and from the northeast sectors (Fig. 3), and are channeled 
into these two main directions by the general atmospheric 
circulation of the region and by the topography formed by 
the Alps and the island of Corsica. Over the period from 
February 1999 to July 2003, only five days were recorded 
with a wind speed above 40kts, and 34 and 100 days with 
wind speeds between 35-40 kts and 30-35 kts, respectively. 
These high wind speeds, and the associated large swells are 
concentrated in the late fall to early spring time periods of 
November to March, respectively. 


North 



South 


Fig. 3. A so-called wind rose for the BOUSSOLE 
site showing the dominant wind directions (the data 
are from the Meteo buoy). 

The range of environmental parameters is also illus- 
trated in Fig. 4, where physical conditions are displayed 
as a four- year record of the wind speed, significant wave 
height, and sea surface temperature (SST). These data are 
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Fig. 5. A time series of the total chlorophyll a (TChl a) concentration at the BOUSSOLE site. Low to 
high concentrations are represented by light to dark colorations in the gray-scale palette, respectively. The 
calendar years for the time series are given along the top of the plot. 


obtained from a weather mooring, which is approximately 
two nautical miles away from the BOUSSOLE site. The 
weather mooring is instrumented and maintained by Meteo 
France, the French weather forecasting agency. 

The minimum SST is about 12.7°C (associated with a 
salinity of 38.4 psu), which is a constant value reached in 
winter when the water mass is fully mixed down to the bot- 
tom. This deep mixing contributes to the formation of the 
bottom waters of the western Mediterranean Sea. These 
data are not further discussed, as they are provided just 
for illustration of the typical annual cycle of the physical 
conditions at the BOUSSOLE site. 

3.1.3 Trophic Conditions 

Oligotrophic conditions prevail at the BOUSSOLE site 
during the summer with chlorophyll a concentrations be- 
low 0.1 mgm -3 , with minima around 0.05 mgm -3 . Higher 
concentrations, up to about 5 mg nr -3 , occur during the 
early spring bloom (February to March or April) . Intermit- 
tent enhancements occur in winter when sunny weather oc- 
casionally, and temporarily, stabilizes the nutrient-rich wa- 
ters. Moderate concentrations, between 0.1-0.2mgm -3 , 
characterize most of the other periods of the year. This is 
illustrated in Fig. 5 and later on in this report (Sect. 8). 

3.1.4 Other Activities in the Area 

The BOUSSOLE site (Fig. 1) is located in an area that 
has been defined some years ago as an area dedicated to 


scientific work, and where is located another site where 
monthly cruises are taking place in the frame of the DY- 
FAMED observation service started in 1991 (Marty 2002f). 
This service collects core data that are made publicly avail- 
able to the entire scientific community and includes CTD 
casts, HPLC pigments, nutrients, oxygen, dissolved or- 
ganic carbon (DOC), and short-time 14 C incubation pri- 
mary production. Other less permanent activities are car- 
ried out occasionally around this site [see the Deep-Sea 
Research special issue 49(11), 2002, Studies at the DY- 
FAMED French Joint Global Ocean Flux Study (JGOFS) 
time-series station, N.W. Mediterranean Sea]. 

3.2 Coastal Transects 

The transect from the coast (actually from the Nice 
harbor) to the BOUSSOLE site is characterized by a per- 
manent front between the waters of the Ligurian Current 
and the more offshore waters. This front moves from near 
the coast (about lOnrni) in spring and summer to near the 
BOUSSOLE site in winter. 

In addition to classic hydrography (CTD casts), which 
have already been performed in the past along the Nice- 
to-BOUSSOLE transect (not systematically, however), the 
BOUSSOLE activity is documenting the IOPs (total ab- 
sorption and attenuation, backscattering at three wave- 
lengths, and CDOM fluorescence) at six fixed stations. 

f Additional information is available at the following Web site: 

http : / /www . obs-vlf r . f r/ sodyf /home . htm. 
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Table 1 . A summary of the types of data collected in support of the BOUSSOLE activity at the coastal, buoy, and ship 
sites (the latter is in very close proximity to the buoy). The above- and in- water references indicate how the oceanic 
data are collected (for example, above- or in- water radiometers). The monthly service visits to the buoy, although 
systematically scheduled, only occur if weather, site conditions (e.g., wave height), and instrument functionality permit. 
The data acquisition rates for the coastal and buoy continuous sampling are as follows, respectively: a) variable, down 
to 10 min; and b) 15 min (lrnin at 6 Hz every 15 min from sunrise to sunset) plus 1 min at 6 Hz each hour during 
nightime. 


Sampling 
(and Site) 

Data 

Type 

Instrument f 
or Model 

Physical 

Measurements 

A 

[nm] 

Data 

Products 

Continuous 

Atmospheric 

CE-318 

Sky radiances, and attenuation of 

A 6 

Aerosol optical thick- 

(Coastal) 



the direct solar beam 


ness (r a ) and types 

Continuous 

Above- Water 

OCI-200 

Downward irradiance 

A7 

E s 

(Buoy) 

In- Water 

OCI-200 

Downward irradiance at 4 and 9 m 

A7 

E d 


In- Water 

OCI-200 

Upward irradiance at 4 and 9 m 

A7 

E u 


In- Water 

OCR-200 

Upwelling radiance at 4 and 9 m 

A7 

L u and Lw 


In- Water 

MINITracka 

Fluorescence at 4 and 9 m 


F, [TChl a] proxy 


In- Water 

C-Star 

Beam attenuation at 9 m 

660 

c 


In- Water 

SBE 37 SI 

Hydrography at 9 m 


T, S, and P 


In- Water 

Hydroscat-II 

Backscattering at 9 m 

A2 

b b 


In- Water 

EZ-III 

Tilt and heading at 9 m 


Tilt and heading]: 

Monthly 

Above- Water 

SIMBADA 

Total water radiance and atten- 

An 

L\y and r a 

(Ship) 



uation of the direct solar beam 




Above- Water 

OCMOOO 

Downward irradiance 

A13 

E s 


In- Water 

oci-iooo 

Downward irradiance 0-150 m 

A13 

E d 


In- Water 

OCI-IOOO 

Upward irradiance 0-150 m 

A13 

E u and L\y 


In- Water 

SBE 911+ 

Hydrography 0-400 m 


T, S, P, F, and [0 2 ] 


In- Water 

Rosette§ 

Discrete seawater samples 0-400 m 


Ci , cip , and 


In- Water 

AC-9+ 

Absorption, attenuation 0-400 m 

A9 

a, c, and b = c — a 


In- Water 

ECO-BB3 

Backscattering coefficient 0-400 m 

A3 

b b 


In- Water 

C-STAR 

Attenuation coefficient 0-400 m 

Al 

c 


In- Water 

WCF1 

CDOM fluorescence 0-400 m 


CDOM description 


Notes: 


f Instrument names, as well as the corresponding manufacturers, are cited here, and in other sections of this report, because 
they have been identified as appropriate for the particular application. This identification is not a recommendation for the 
use of any of these instruments; nor does it imply that these instruments are the best suited for the application or for other 
similar projects. In addition, the possibility of using any other instrumentation than those cited here is retained without 
notice. 

I The buoy heading is used to determine the orientation of the in-water sensors, which are mounted on fixed arms, with 
respect to the Sun. 

§ SBE 12-bottle rosette for laboratory analyses of phytoplankton pigment concentrations (Ci) from HPLC; plus particles and 
phytoplankton absorption coefficients (a p and a respectively). 

WETLabs CDOM Fluorometer. 

Ai 660 nm. 

A 2 442 and 560 nm. 

A 3 470, 530, and 660 nm. 

Ag 440, 675, 870 (plus 870 polarized), 940, and 1,020 nm. 

A 7 412, 443, 490, 510, 560, 665, and 683 nm (a spare set of sensors have 443, 490, 510, 555, 560, 665, and 683 nm). 

Ag 412, 443, 488, 510, 532, 555, 630, 676, and 715 nm. 

An 350, 380, 412, 443, 490, 510, 565, 620, 670, 750, and 870 nm. 

A13 412, 443, 455, 490, 510, 530, 560, 620, 665, 683, 705, 780, and 865 nm before July 2003, and 380, 412, 443, 455, 490, 510, 

530, 560, 620, 665, 683, 705, and 780 nm thereafter. 
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The stations are approximately 5nmi away from one an- 
other, and the survey work started in May of 2003. Docu- 
menting these properties at this series of stations is also a 
way to put the observations made at the BOUSSOLE site 
into a broader context. 

3.3 The Coastal AERONET Site 

Aerosol amounts and types are information of prime 
importance in the process of vicariously calibrating ocean 
color sensors (Sect. 10), but they cannot be obtained easily 
at the mooring site. This is the main reason for having 
installed an automatic scanning sun and sky photometer at 
Cap Ferrat , which is providing a continuous record (cloud 
permitting) of these parameters thanks to the logistics and 
data transmission capabilities of the AERONET. 

The parallel analysis of the sun photometer and satel- 
lite observations are expected to indicate the circumstances 
wherein the data collected at the coast are representative of 
the conditions offshore. In support of this study, aerosol — 
and when the sky is clear — optical thickness measurements 
are also performed en route when steaming from Nice to 
the BOUSSOLE site, in order to get some idea of the spa- 
tial heterogeneity of aerosol properties. 

4. DATA COLLECTION 

The three subsections below briefly summarize the set 
of measurements that are performed offshore, either from 
the ship or from the buoy, and at the coastal site. They 
are summarized in Table 1. Details on the data acquisition 
protocols will be provided in Sect. 6 and data processing 
in Sect. 7. 

4.1 Ship Measurements 

The set of parameters directly derived from the mea- 
surements made by the above-water SeaWiFS Multichan- 
nel Surface Reference (SMSR) and the in-water SPMR to 
200 m depth are as follows (A 13 denotes the 13-channel 
spectral band set): a) surface irradiance, E s ( A 13 ), and b) 
upward and downward irradiance, E u ( A 13 ) and E d ( A 13 ), 
respectively. The wavelengths for these measurements are 
380, 411, 443, 456, 491, 510, 532, 560, 620, 665, 683, 705, 
779, and 865 nrn 

Profiles with the CTD rosette to a depth of 400 nr are 
used to collect the following parameters: chlorophyll fluo- 
rescence, CDOM fluorescence, attenuation and absorption 
coefficients (412, 440, 488, 510, 532, 555, 630, 676, and 
715 nrn), (440, 532, and 650), temperature, conductiv- 
ity and derived salinity, oxygen concentration, and surface 
PAR. Seawater samples are collected and filtered for sub- 
sequent determination of the phytoplankton pigments and 
particulate absorption. 

The SIMBADA above-water radiometer collects water- 
leaving radiance and aerosol optical thickness at 350, 380, 
412, 443, 490, 510, 565, 620, 670, 750, and 870 nrn. 


4.2 Buoy Measurements 

The set of parameters directly derived from the buoy 
above- and in-water measurements are as follows (A 7 de- 
notes the 7-channel spectral band set for the AOP sensors): 
a) E s ( A 7 ); b) E d ( A 7 ), E u ( A 7 ), L u ( A 7 ), c(660), and chloro- 
phyll fluorescence at 4m; and c) Ed{ A 7 ), E u ( A 7 ), L u { A 7 ), 
chlorophyll fluorescence, c(660), 6^(443), and 6^(560), con- 
ductivity, temperature, pressure, and two-axis tilt as well 
as compass heading at 9 m. The wavelengths for the AOP 
measurements are 412, 443, 490, 510, 560, 665, and 683 nrn 
(a spare set of sensors has the 412 nrn channel replaced with 
555 nrn) . 

From these measurements, various AOPs or IOPs are 
derived: the diffuse attenuation coefficients for upward and 
downward irradiance, K u and K d \ the attenuation coeffi- 
cient for upwelling radiance, Kp\ the diffuse reflectance 
just below the sea surface, R ; the nadir Q-factor, E u /L u ; 
and the attenuation and backscattering coefficients, c and 
bb. The absorption coefficient, a, is tentatively derived 
through inversion of the AOPs (using, for example, K d 
and R). Here it is assumed that K d and R can be ac- 
curately derived from the two measurement depths of the 
buoy, e.g., Gordon and Boynton (1997), Leathers and Mc- 
Cormick (1997), and Barnard et al. (1999). Table 2 shows 
how the data acquired on the system fits to the band sets 
of several ocean color sensors. 

Table 2. The MERIS wavelengths, as acquired on 
the buoy (leftmost column), and the spectral shift 
of other satellite bands with respect to those of 
MERIS. The second row gives the channel band- 
widths. 


MERIS 
10 nrn 

POLDERS SeaWiFS MODIS 
20 nrn 20 nrn 10 nrn 

412 nrn 

443 

490 

510 

560 

665 

683 

0 0 

0 0 0 

0 0-2 
0 

+5 -5 -9 

+5 +5 +2 

0 


The data are derived from the IOCCG (1998) report. 


4.3 Coastal Measurements 

The parameters collected at the coastal site are the sun 
and sky radiances, using standard Liters at 440, 675, 870, 
940, and 1,020 nrn plus a polarized Liter for 870 nrn, which 
are used to derive total column water vapor, ozone, and 
aerosol properties, including Angstrom coefficients, optical 
thickness, and polarization ratios! • 


f These derived parameters are available at two Worldwide 
Web addresses: http://www-loa.univ-lillel.fr/photons 
/DATA/Villefranche and http://aeronet.gsfc.nasa.gov. 
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5. THE BOUSSOLE BUOY 

The description of the BOUSSOLE buoy is presented 
more or less in a chronological order, starting with the 
design, and then continuing with the successive versions 
and the various tests that were performed. 

5.1 Design Concept 

Platforms developed for oceanographic purposes are 
rarely adapted to the unique problems of deploying ra- 
diometers. Accurately measuring the in- water light held 
is difficult, because the instruments themselves and, usu- 
ally more significantly, the platform onto which they are 
installed, introduce perturbations (shading in particular). 
Other difficulties originate from the need to keep the in- 
struments oriented as stably as possible, either because a 
nadir radiance or a planar irradiance is being measured. 
The actual measurement depth can also be difficult to ac- 
curately assess, because rapid vertical displacements of the 
instruments sometimes occur, which prevent any precise 
estimation of pressure (and, thus, depth). Considering the 
above difficulties (among others), a new type of platform 
was developed and specifically dedicated to making high- 
quality radiometric measurements at sea. The platform 
minimizes shading effects by reducing the cross-sectional 
area of structural components, and its wave-interaction 
characteristics ensure the stability of the instruments. 

The design constraints for the new platform were: 

1. Measure E u , Ed , and L u (nadir) at two depths, plus 
E s at the surface; 

2. Minimize the shading of the instruments; 

3. Maximize the stability of the instruments; and 

4. Permit deployment at a site with a water depth of 
2,440 m, and swells up to 8 m (but small horizontal 
currents) . 

The basic design principle for the buoy is that of a re- 
versed pendulum , with Archimedes thrust replacing grav- 
ity. A large sphere (with a diameter of about 1.8 m) is sta- 
bilized at a depth out of the effect of most surface waves, 
and connected at the end of a long cable anchored on 
the sea bottom. This sphere creates the main buoyancy 
of the system. A rigid, tubular, structure is fixed above 
the sphere, which hosts the instrumentation on horizontal 
arms (at 4 and 9m depths). The resulting approximately 
three tons of thrust ensures the stability of the system, 
which is subjected to very limited forces from the so-called 
transparent-to-swell superstructure. 

The taut-cable mooring is significantly different from 
what is usually referred to as spar buoys. With the BOUS- 
SOLE design, there is no large body at the surface gener- 
ating shade, and the stability of the instruments is ensured 
even for quite large swells. The small vertical excursions of 
the near-surface structure means the water level above the 
instruments is accurately known at all times. A schematic 
depiction of the buoy design is presented in Fig. 6. 



Solar 

Panel 



Fig. 6. A schematic of the principle BOUSSOLE 
components distinguished according to the upper 
and lower superstructure portions (total height is 
25 m) . The in- water radiometers are mounted at the 
end of the sensor arms which are located well above 
the buoyancy sphere (the largest perturbation to 
the light held; the underwater structure is entirely 
black, however). The above- water solar reference is 
sited above the solar panels, which are located 4 m 
above the nominal seawater surface. 
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Fig. 7. The complete mooring line on the deck of the Castor-02, before deployment. The beginning of 
operations started 20 July 2000, 0700 Greenwich Mean Time (GMT). 


Theoretical calculations were made on the response of 
the proposed design to the prevailing environmental condi- 
tions, by specifying an initial and preliminary design and 
material for the construction, as well as a swell height 
of 5 m and a 7 s wave period (typical of the deployment 
site). The natural response frequencies of the whole moor- 
ing were determined along with the period and amplitude 
of the oscillations and displacements due to swell and cur- 
rents. These calculations were extremely encouraging in 
terms of the expected tilt and oscillation parameters, so 
the decision was made to first build a reduced-scale model 
in order to perform more realistic tests. 

5.2 Reduced-Scale Model 

A one-tenth scale model was developed in order to 
verify the theoretical predictions, and it was tested in a 
24 x 16x3 nr 3 engineering pool with a 10 m shaft, by ap- 
plying to it several types of swells and currents, either 
monochromatic or random swells, up to 5 m real scale. 
Results of these tests completely confirmed the theoreti- 
cal predictions that were previously carried out, in terms 
of horizontal and vertical displacements (the latter are ex- 
tremely low) , as well as in terms of angular deviations from 
the vertical. For instance, the mean tilt of the buoy is ap- 
proximately 4° (with about ±4° of pitching), for a 4.6m 
swell of period 5.2 s, i.e., at the limit before wave breaking 
occurs. 

The reduced-scale model tests also confirmed no hid- 
den defects were present — which are usually undiscernible 


with calculations alone — so there was no reason to believe 
there was an unknown element that might compromise the 
feasibility of the overall system. The tests also showed a 
significant sensitivity to currents, but the deployment site 
is typified by rather low current velocities. The tests were 
considered fully conclusive and led to the construction of 
the first full-scale version 1 of the buoy. 

5.3 Full-Scale Version 1 

The version 1 buoy was built during the spring of 2000. 
To minimize weight, aluminum was used as the principal 
construction material, and the buoy was manufactured in 
two large parts. The lower section, which extends from 
20 m to 9 m below the surface, consists of the buoyancy 
sphere and a central tubular support structure; the upper 
section, which goes from 9 m below the sea surface to 4.5 m 
above the surface, hosts all of the above- and in- water sci- 
entific instrumentation either within the extension of the 
central tubular structure or at the ends of sensor arms. 

The so-called beta (version 1) buoy was deployed dur- 
ing three months (20 July to 20 October 2000), which 
exposed the system to a variety of meteorological situa- 
tions. The buoy was equipped with two inclinometers, a 
pressure sensor, an ARGOS beacon, and a flashing light. 
The goal of this deployment was to evaluate the design 
concepts of the buoy so as to identify possible problems 
requiring modifications. Figure 7 shows the full mooring 
on the deck of the ship before being deployed. A full day 
was necessary for the deployment operations, which were 
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significantly different than what is usually performed for 
classical moorings because it is a taut mooring. The de- 
tails of these operations are not given here. The results of 
the initial deployment are shown in Figs. 8 and 9. 



Fig. 8. The percentage of the tilt measurements 
that were less than 5° , as computed over successive 
10-day periods of the qualification deployment. 

Ignoring details beyond the scope of the initial deploy- 
ment, the following points are relevant: 

■ From mid-September to the end of the deployment, 
the conditions were quite rough (mean swells up to 
4 m with at least 5 m peaks) . 

■ Until mid-September, the periods where the buoy 
is inclined at angles larger than 5° is small, which 
means that the beta version was already satisfying 
the requirements for meteorological conditions from 
calm to moderate. 

■ Oscillations were observed, which exactly match the 
inertial period of the ocean at the latitude of the 
site, and are most likely due to the frequency re- 
sponse of the water masses after the wind stress is 
relaxed. 

In summary, the results of the qualification deployment 
showed some lack of righting torque, with discernible con- 
sequences on the behavior of the buoy when the swell be- 
comes greater than about 2 m. Although not dramatic, 
this behavior is unsatisfactory, so slight modifications in 
the design and construction of the buoy were introduced 
in order to improve the percentage of time for which the 
requirements — in terms of inclination — are satisfied. 

5.4 Version 2 Modifications 

For the version 2 design of the buoy, the buoyancy 
sphere was placed at 17 instead of 18 m. Consequently, the 
distance from the center of the sphere to the termination 


point of the buoy and the mooring cable was 1 m greater 
(and the righting torque was, therefore, much larger). The 
materials used to build the upper part of the buoy (the one 
equipped with arms and the instrumentation) were also 
changed from aluminum to carbon composite. A specific 
study was conducted to verify the composite material was 
a satisfactory substitute, and this was definitely the case. 
This study was also used for dimensioning the material 
used during construction. 

The impact of these changes on the upper part of the 
system were rigidity was improved, weight was reduced by 
approximately two-thirds, corrosion risks were decreased, 
the drag coefficient was reduced by the use of superstruc- 
ture tubes with a smaller diameter, and the elastic limit 
was larger than it was in comparison to the forces imposed 
by breaking waves. In addition, a net gain of 60% in the 
righting torque was obtained for the full buoy. 

5.5 May 2002 Deployment 

The modified version of the buoy was deployed onsite 
with all its instrumentation on 16 May 2002. Some images 
of this deployment are shown in Fig. 10. Soon after this 
deployment, an atypical behavior of the buoy was identi- 
fied, in particular, an abnormal lowering of the buoy — the 
exposed superstructure was too low in the water. When 
a clear analysis of what was happening was possible, and 
when the weather permitted safe at-sea operations, a ship 
was sent onsite. 

Unfortunately the ship arrived too late, and the entire 
system was lost. The cause of the sinking was a leak in 
the main flotation sphere. The exact reason for the leak 
was not discovered, but most likely it was caused by a 
manufacturing defect. 

5.6 Version 3 Modifications 

After the May 2002 accident, two independent engi- 
neering studies were enacted, the objective of which was 
to improve the confidence in the design of the buoy and 
to provide advice for the construction of a replacement. 
These two studies were performed by the Institut Frangais 
de Recherce pour V Exploitation de la Mer (IFREMERf) 
and the Marine Technology Research Institute): (MARIN- 
TEK). The results of these inquiries confirmed the feasibil- 
ity of the transparent-to- swell concept, but recommended 
the use of either foam or steel for the construction of the 


f Translated as the French Research Institute for Exploitation 
of the Sea and headquartered in Issy-les-Moulineaux, France. 
Additional details are available at the following Web address: 

http : //www. if remer . fr. 

i Headquartered in Trondheim, Norway, MARINTEK devel- 
ops and verifies technological solutions for the shipping and 
maritime equipment industries and for offshore petroleum 
production. Further information is available at the following 
Web address: http://www.marintek.sintef.no 
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Fig. 9. The results of the qualification deployment (20 July to 20 October 2000), showing the time series of 
a) the depth of the buoy CTD (nominally at 8.3 m when the buoy is motionless and the sea surface perfectly 
flat), b) the buoy tilt, c) the swell height (the one-third most significant waves), d) the swell period, and e) 
the wind speed. 
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Fig. 10a. The upper superstructure of the buoy (Fig. 6) prior to it being fastened to the lower superstructure, 
which is done under water by divers. The optical sensors are visible at the end of the sensor arms in the left 
half of the picture. The two oblique arms are only used during deployment, to secure the buoy arms; they 
are removed after the buoy is deployed. 



Fig. 10b. The buoy completely deployed with the top-most portion of the structure rising approximately 
4m above the water surface (as measured from the bottom of the solar panels). 
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buoyancy sphere. The second option was selected, and a 
special design was prepared, with a segmentation of the 
total volume into seven pressure-resistant and watertight 
compartments. This design minimizes the risk of rapidly 
loosing a large part of the buoyancy in case a leak occurs. 
The upper superstructure is still built with aluminum. 

In addition, a strain gauge was installed at the base 
of the buoy, just at the beginning of the mooring cable. 
Because the tension in the cable that is recorded by the 
gauge is, by definition, equal to the buoyancy of the entire 
mooring, any change in buoyancy can be detected as a 
change in tension. 

5.7 Fall 2003 Deployment 

The new version 3 buoy and mooring were deployed 6 
September 2003, and the upper instrumented part was re- 
covered three months later (6 December). The behavior of 
the buoy was as expected, and the data collection was con- 
tinuously operating during the three-month period, with 
the exception of a) the transmissometer, which flooded af- 
ter about lOdays (because of a faulty o-ring seal of the 
pressure housing); and b) the strain gauge (which was sub- 
sequently replaced with a new sensor). 

Several storms occurred during the deployment time 
period. In particular, one had a wind speed of 50kts, 
gusts up to 65kts, and maximum wave heights of about 
8 m (such a severe sea state never occurred within the past 
five years). In such conditions, the strong wind-driven sur- 
face current pushes the buoy, which starts to tilt and si- 
multaneously sinks, because the Kevlar mooring cable can- 
not stretch. After several hours of a strong wind, the buoy 
simply becomes totally submerged, with the top-most por- 
tion of the buoy reaching maximum depths of about 5 m. 
This situation is not risky for the instrumentation, which 
is designed to go even deeper, and is actually an effective 
mechanism for protecting the buoy from breaking waves, 
which are the only potential source of severe damage. Af- 
ter the wind stops, the buoy rapidly rights itself and goes 
back to its nominal position. 

The upper superstructure severely suffered from corro- 
sion, the cause of which was the use of improper sacrificial 
anodes. The faulty anodes were replaced by anodes specifi- 
cally manufactured for aluminum structures, and corrosion 
was no longer a problem. The upper superstructure used 
for this deployment was not considered reusable, so a new 
one was built in preparation of the next deployment. 

5.8 March 2004 Deployment 

The upper superstructure of the buoy, which hosts the 
instrumentation, was reinstalled on 4 March 2004; this 
should have been done on 4 February, but a last-minute 
problem with the buoy computer postponed the deploy- 
ment, and eventually led to a one-month delay. The de- 
ployment lasted until 23 July 2004, when the rotation of 


the upper superstructure took place. Indeed, during the 
first five months of 2004, a second buoy superstructure 
was built and a second suite of instrumentation was ac- 
quired, in order to allow a bi-yearly rotation of the full 
equipment. The rotation was performed by combining a 
ship and divers on site and transportation of buoys from 
the site to the coast (and vice versa) by helicopter. 

5.9 July 2004 to June 2005 Deployment 

This deployment lasted 10 months (end of July 2004 
to beginning of June 2005). Only a few difficulties arose, 
mostly due to stops of the acquisition system, requiring 
replacement by divers of the main buoy computer. 

5.10 Full Rotation (June 2005) 

During the end of 2004 and the beginning of 2005, a sec- 
ond lower superstructure and a second mooring line were 
built, in order to prepare for a complete rotation of the 
entire mooring line and buoy. This operation was pri- 
marily motivated by the unknown longevity of the 2 km 
long Kevlar mooring cable. A 21-month deployment was 
deemed sufficient for a first check to be performed. An- 
other motivation was to check the integrity of the main 
buoyancy sphere. The operations took place from 2-3 June 
(recovery of the mooring and buoy) and 15-16 June (de- 
ployment of the new system) . The system has been at sea 
since then. 

6. INSTRUMENTATION 

The scientific instrumentation are categorized accord- 
ing to the deployment site: buoy (which includes the ac- 
quisition system), ship, and coastal. The description of the 
instruments are separated into sections explaining how the 
radiometers are characterized and how bio-fouling is min- 
imized on the buoy sensors. 

6.1 Buoy Sensor and Acquisition System 

In this section, brief descriptions of the buoy sensors 
and the acquisition system are provided. 

6.1.1 Sensors and Data Collection Scenario 

Instruments on the buoy, simultaneously collecting data 
in a continuous way, include the following: 

• Satlantic ocean color 200-series radiometers mea- 
suring Ed, E u , and L u (nadir) at two depths (4 and 
9 m); 

• A Satlantic Multichannel Visible Detector System 
(MVDS) 200-series radiometer measuring E s (at 4.5 
above the sea surface); 

• An Advanced Orientation Systems, Inc. (Linden, 
New Jersey) two-axis tilt and compass sensor at 9 m 
(EZ-III); 
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• A Sea-Bird Electronics, Inc. (Bellevue, Washing- 
ton) CTD at 9 nr measuring conductivity, tempera- 
ture, and pressure. 

• Chelsea (Surrey, United Kingdom) MINItracka Mk 
II fiuorometers at 4 and 9 nr for a proxy for the 
TChl a concentration, [TChl a] . 

• WETLabs C-star transnrissonreters at 4 and 9 nr for 
the beam attenuation coefficient at 660 nnr, a proxy 
of the particle load. 

• A Hydro-Optics, Biology, and Instrumentation 
Laboratories (HOBI Labs), Inc. (Tucson, Arizona) 
Hydroscat-II backscattering meter at 9 nr measur- 
ing a proxy to 6&( A) at two wavelengths (442 and 
560 nnr) . 

These data are collected every 15nrin during daylight, and 
at hourly intervals at night. Each data acquisition se- 
quence lasts one minute. 

6.1.2 The Integrated Acquisition System 

The central recording device of the buoy is the Data 
Aquisition and Control Network (DACNet) unit, which 
houses the primary computer: a PC104 Cool RoadRunner 

11 200 MHz 6 x86 with 32 MB of random access memory 
(RAM) and the serial data acquisition equipment (record- 
ing to a 1 GB disk drive) . These components are used in 
the collection, storage, and downloading of the data ob- 
tained from the instruments. The majority of the time the 
system is inactive, and the DACNet unit is unpowered. A 
small internal microcontroller is responsible for the power 
supervision and control, and a precision clock. 

Optimization of the power budget is an important el- 
ement, because of the restricted power available from the 

12 VDC 105 Ah underwater battery, which is recharged us- 
ing solar panels. The power managment functions can be 
summarized as follows: 

1. An accurate real-time clock; 

2. An alarm clock for powering up the system for the 
programmed schedule; 

3. A watchdog timer that protects against draining the 
battery in the event of a computer malfunction; and 

4. A power-fail shutdown to prevent computer startup 
if the battery power is too low. 

The DACNet operating system runs on Red Hat Linux 

6.2 kernel v2.2.19, which accommodates custom-made Java 
software: the Satlantic Telemetry Acquisition Manager 
(STAM) and the Satlantic Node Manager. These pro- 
grams provide the operational, configuration, and interface 
requirements for the buoy to run autonomously and enable 
communication of data and hies with the user. 

A third piece of software, Satlantic Base Manager, runs 
on a personal computer (PC) equipped with the Windows 
NT or Windows 98 operating system and provides the op- 
erator with administrative tools — each with a graphical 


user interface (GUI). These tools enable the user to con- 
figure the Node Manager, transfer data and hies, update 
the internal clock of the acquisition node, view data from 
specihed instruments in real time, and completely shut- 
down the acquisition node for maintenance!- 

6.1.3 Data Communications 

The complete data stream can be downloaded through 
a wireless ethernet link established between the buoy and 
the ship (the same equipment is used aboard the buoy and 
the ship). This link functions with the ship positioned 
a few hundred meters from the buoy. The connection is 
driven by the buoy, with a wakeup sequence each hour 
of the day. This hierarchy was adopted, rather than es- 
tablishing the communication from the ship, in order to 
minimize the time during which the buoy communication 
hardware is awake and consuming energy. 

Data retrieval from the buoy is performed from the 
DACNet via a Cisco Aironet 340 series wireless bridge 
manufactured by Cisco Systems (San Jose, California) and 
is stored directly on a PC. The data is retrieved in a bi- 
nary format and consists of one daily log hie for each 
of the following seven groups: a) the instruments con- 
nected to the DATA-100 at 4 m (radiometers and fiuorom- 
eter); b) the instruments connected to the data acquisi- 
tion module (DATA-100) at 9 m (radiometers, huorometer, 
transmissometer) ; c) the CTD; d) the Hydroscat-II; e) the 
strain gauge; f) the MVDS; and g) the two-axis tilt sensor 
and compass. The daily hie is closed once communication 
through the Cisco system occurs. In this case, a new hie is 
created starting from the next measurement taken by the 
DACNet system until midnight or until a subsequent link 
up via the Cisco system is performed. Storage space on 
the hard disk of the DACNet has the capacity to store up 
to approximately three months of data. 

Part of the data stream is transmitted via the ARGOS 
system and is used for surveying the functioning of the 
system; the sample data include the tilt and depth of the 
buoy, the strain of the mooring cable, the battery voltage, 
the disk space, the spectrum of the above- water irradiance, 
and instrument health parameters, which are indicating 
whether or not instruments and the acquisition system are 
functioning nominally. 

6.2 Shipboard Sensors 

In this section, brief descriptions are provided of the 
instruments deployed from the ship to help characterize the 
BOUSSOLE site during the monthly maintenance cruises. 
The objective here is to have quality-assured data for the 
most important hydrographic, optical (AOPs and IOPs), 
and biogeochemical parameters. 


f More details of this software can be found in Satlantic’s man- 
ual DACNet Software Overview — Villefranche Remote Opti- 
cal Mooring. 


17 



BOUSSOLE: A Joint CNRS-INSU, ESA, CNES, and NASA Ocean Color Calibration and Validation Activity 


6.2.1 CTD and Water Sampling 

The CTD/rosette package is based on a Sea-Bird Elec- 
tronics (SBE) 911 plus, equipped with the following sen- 
sors: 

■ Pressure (Digiquartz Paroscientihc, Inc.); 

■ Temperature (SBE 3); 

■ Conductivity (SBE 4); 

■ Dissolved oxygen (SBE 13 until January 2003, and 
then an SBE 43); and 

■ Altimeter (Datasonics PSA 900). 

The sensors are mounted on a SBE 32 rosette with 11 12 L 
Niskin bottles (the 12th being replaced by an AC-9+ as 
discussed below). 

The pressure sensor is not regularly calibrated (it was 
last calibrated in 1995); it is known from past experience 
with this type of sensor that the uncertainty is on the order 
of 0.5dbar. It is acknowledged that this type of sensor 
only experiences a drift of the zero point, which is easy 
to control, and to correct for, when the sensor is at the 
surface. 

The temperature sensor is sent once a year for calibra- 
tion to SBE, which has the agreement from the National 
Institute of Standards and Technology (NIST) for this type 
of work. Figure 11 provides an example of the calibrations 
performed in January 2002 and January 2003 (was about 
the same between January of 2001 and January of 2002). 
The drift is on the order of 0.0025° C over one year, and it 
is linear. Any temperature value computed between two of 
these calibrations is, therefore, provided with a ±0.0025°C 
uncertainty. 
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Fig. 11. The calibration drift of the temperature 
sensor mounted on the CTD. 

The conductivity sensor is sent once a year to SBE for 
calibration. Additional calibration against values deter- 
mined from water samples and subsequent laboratory anal- 
yses (with a laboratory salinometer) are not performed. 
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Figure 12 provides an example of the calibration performed 
in 2001. The change in the calibration slope is caused by 
biological fouling that grow inside the sensor, modifying 
its geometry. The influence of this fouling is larger as the 
conductivity increases. Typical values of the BOUSSOLE 
site are on the order of 5Sm _1 . It is, therefore, estimated 
that the sensor drift over one year is on the order of 0.0005 
in conductivity units, which corresponds to approximately 
0.005 psu. The drift was a bit higher in 2002, around 0.001 
in conductivity units, which is about 0.01 psu. 
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Fig. 12. The calibration drift of the conductivity 

sensor mounted on the CTD. 

From July 2001 to December 2002, a SBE 13 sensor 
with a Beckman membrane was used. This sensor showed 
a large drift. Starting in January of 2003, a new sensor 
was used (SBE 43), which is provided with a nominal 5% 
uncertainty and a less than 2% drift. No Winkler-based 
calibrations are performed on these sensors. 

At each cast, a first round trip down to 50 m is per- 
formed with the CTD, in order to rinse all instruments 
and equilibrate them in temperature. Then the package is 
lowered down to 400 nr at a descent rate of about 0.5 nrs -1 . 
The data are collected at a 24 Hz frequency, and then pro- 
cessed through the SeaSoft software until binned values are 
produced at a 1 nr resolution. Each data hie has an asso- 
ciated header hie containing the necessary ancillary infor- 
mation (latitude, longitude, date and time of acquisition, 
meteorological information, etc.). 

Some of the sampling depths are hxed (i.e., 5 and 10 nr 
for surface analyses) and others are determined at the end 
of the rosette descent, based on the phytoplankton flu- 
orescence prohle, in order to best represent any specihc 
features of the vertical prohle. No sampling is done deeper 
than 200 nr. Bottles are closed during the ascent of the 
rosette. When back on deck, 2.8 L plastic bottles are used 
to collect water that is used for subsequent hltration. The 
hlters are kept in Petri slides, wrapped in aluminum and 
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stored first in liquid nitrogen on board and then in a — 80° C 
freezer in the laboratory. The material collected on filters 
are used for HPLC analyses of the phytoplankton pigments 
(Sect. 7.1) and filter pad absorption (Sect. 7.2). For the 
surface samples, triplicate hltrations are performed for the 
HPLC analyses. 

6.2.2 Fluorometer 

A Chelsea Mk III Acquatracka fluorometer is installed 
on one of the analog channels of the CTD SBE 911. This 
instrument is equipped with an excitation filter at 430 nrn 
and an emission filter at 680 nrn. The data are provided 
both as voltages and chlorophyll concentrations, the latter 
being computed with the initial calibration coefficients pro- 
vided by the manufacturer. These concentrations are just 
indicators. Another calibration of the fluorometer signal 
can be performed against the HPLC discrete chlorophyll 
concentration values. 

6.2.3 AC-9 

One of the CTD rosette bottles is replaced by a 25 cm 
path AC-9+ sensor. This instrument measures the absorp- 
tion and attenuation coefficients of the water at nine wave- 
lengths (Table 1). In order to visualize and collect the data 
in real time, the A09+ is linked to the surface through a 
19,200 baud serial connection, which also provides power to 
the instrument. Although the CTD and A09+ can be in- 
dependently controlled, only one computer is used in order 
to provide the same time frame and depth measurements 
to both the hydrographic and optics data. The connec- 
tion only functions in the computer-to-AC~9+ direction, 
although theoretically two-way communications are possi- 
ble. The choice of one-way communications is driven by 
the different transmission speeds (300 baud from the com- 
puter to the SBE 911 and 19,200 baud from the SBE 911 
to the AC-9+). 

An SBE 25 pump is installed at the entrance of the 
AC~9+ water circuitry, and it is coupled with the pump 
that provides water to the CTD conductivity sensor. Both 
only function when submerged, which is detected when 
a certain conductivity threshold is reached. The AC-9+ 
is programmed to start the following acquisition sequence 
immediately after being powered by the CTD sensor: 15 s 
time lag before starting acquisition, then a 60s warm- 
ing period, then a 60s period allowing for cleaning of the 
plumbing components, and then a 35 min data collection 
period. The latter is enough for the round-trip down to 
400 m, including the stops for bottle firing. Note that the 
starting sequence of 135 s is occurring during the initial 
round trip down to 50 m. This procedure ensures a pro- 
gressive and full temperature rise of the electronics and a 
correct water circulation into the instruments. Data col- 
lection is executed using the WETLabs WetView software 
v5.0a (as provided by the manufacturer). 


6.2.4 ECO-BB3 

A WETLabs three-channel backscattering (ECOBB3) 
meter is connected to one of the auxiliary ports on the 
AC~9+ (a serial port with a 19,200 baud rate). This sen- 
sor emits light at 440, 532, and 650 nm, and measures the 
backscattered light at the same three wavelengths. Power 
to the instrument is also provided by the CTD. The data, 
however, cannot be visualized in real time; they are stored 
in the AC-9+ memory, and then collected once the rosette 
is back on deck and the 35 min acquisition sequence of the 
AC~9+ is terminated. 

6.2.5 CDOM WETStar Fluorometer 

A WETLabs WETStar CDOM fluorometer is plugged 
into one of the analog ports of the CTD. This instrument 
has an excitation wavelength at 370 nm and an emission 
filter at 460 nm. It is supposed to provide a qualitative 
measure of the CDOM concentration. It is connected to 
the plumbing used for circulating water through the CTD 
sensors; the water goes from the conductivity sensor to the 
dissolved oxygen sensor, then to the CDOM fluorometer, 
and finally to the pump. The data are collected in paral- 
lel to the CTD data and processed in real time using the 
SeaBird SeaSoft software. The raw signal is provided as 
well as a CDOM concentration computed with the calibra- 
tion coefficients provided by WETLabs. 

6.2.6 SPMR and SMSR 

The Satlantic SPMR free-fall radiometric profiler has 
been used at LOV since 1995, and in many other laborato- 
ries around the world. The LOV version of this instrument 
is equipped with two irradiance sensors, collecting the up- 
ward and downward plane irradiances ( E u and Ed, respec- 
tively) at the following 13 wavelengths: 412, 443, 456, 490, 
510, 532, 560, 620, 665, 683, 705, 779, and 865 nm (the 
latter was replaced by 380 nm after July 2003). 

The SPMR Profiler is made of a long pressure case 
(1.2 m long and 9 cm in diameter) that contains the major- 
ity of the system electronics, while the optical sensors are 
located and separately housed at either end of the cylinder 
and facing in opposite directions (i.e. , up and down in the 
normal descent orientation). One end of the instrument 
has buoyant fins to stabilize the free-fall descent of the 
rocket-shaped device, and the other end has a small an- 
nular lead ballast to establish the righting moment of the 
instrument and to further stabilize the orientation. The 
amount of lead used provides for a fine tuning of the free- 
fall velocity. 

The light sensors used for measuring irradiance (in units 
of milliwatts per square centimeter per nanometer), have 
a black Delrin plate on the end. The plate contains 13 
specially-designed, diffuser-based, cosine collectors. Tilt 
and pressure are recorded at the same frequency as the 
irradiance measurements (6 Hz). 
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The SPMR is accompanied by a deck reference sensor 
(SMSR) which is equipped with the same 13 wavelengths 
as the SPMR and is based on the same electronics. Data 
acquisition is synchronized between the SPMR and the 
SMSR and both operate at the same 6 Hz frequency. 

The absolute calibration (Sect. 6.4) of the SPMR and 
SMSR with respect to NIST-traceable standards is per- 
formed every six months in the Satlantic optics calibra- 
tion laboratory. Temporal stability between calibrations 
is monitored using a portable light source (Hooker and 
Aiken 1998) specifically developed for that purpose by Sat- 
lantic, i.e., the second-generation SeaWiFS Quality Moni- 
tor (SQMTI). Combining these two elements allows a 3% 
maximum uncertainty to be maintained on the calibration 
of the SPMR and SMSR. 

A SPMR profile starts when the instrument reaches 
a distance of 50 m off the stern of the ship (the ship is 
25 nr long and a mark on the cable indicates the 50 m dis- 
tance). The instrument is then released and falls at ap- 
proximately 0.5 ms -1 . The descent is stopped when the 
pressure sensor indicates a depth of about 150 m (except 
in extremely clear waters wherein the profile is performed 
down to 200 nr) . This technique ensures data collection out 
of the ship shadow, and to obtain measurements with tilt 
angles less than 2°. The Sun is usually on the port side 
of the ship, which is not so important because the ship 
shadow is not affecting the measurements. 

6.2.7 SIMBADA 

The SIMBADA instrument is an above- water radiome- 
ter designed and manufactured by the University of Lille, 
France. It is an upgraded version of the SIMBAD (De- 
schamps et al. 2004) above-water radiometer. It measures 
both water-leaving radiance and aerosol optical thickness 
in 11 spectral bands centered at 350, 380, 412, 443, 490, 
510, 565, 620, 670, 750, and 870 nnr by viewing the ocean 
surface (ocean-viewing mode) and the Sun (sun-viewing 
mode) sequentially. The same optics (3° fleld-of-view) , in- 
terference Liters, and detectors are used in both ocean- 
and sun-viewing modes, but different electronic gains are 
used for each mode. The optics are fitted with a vertical 
polarizer, to reduce reflected skylight when the instrument 
is operated in ocean-viewing mode. A small GPS is fixed 
in the front panel of the SIMBADA for automatic acquisi- 
tion of geographic location at the time of data collection. 
Viewing angles are recorded automatically. 

Viewing of the ocean must be made in clear-sky con- 
ditions (3/4 of sky cloudless, and no clouds obscuring the 
Sun), outside the sun glint region (relative angle between 
solar and viewing directions of 45-90°), and at a nadir 
angle of approximately 45°. For those angles, reflected 
skylight is minimized, as well as residual ocean polariza- 
tion effects. The measurements can be made on a vessel 
underway, so there is no need to stop the ship to make 


measurements. To normalize water-leaving radiance, inci- 
dent solar irradiance is not measured, but computed using 
the aerosol optical thickness. 

The operator can select, in addition to ocean- and sun- 
viewing modes, dark current, and calibration modes. Each 
series of measurements lasts 10 s. The frequency of the 
data sampling is approximately 8 Hz. The recorded data 
are stored internally and downloaded using a computer 
serial link at the end of the day or a cruise. The instrument 
is powered by batteries that can be charged using a main 
supply of 110-240 V, 50-60 Hz. About 2h of charging is 
usually enough for one day of measurements. 

The SIMBADA is used from the bow of the ship or 
from the upper superstructure when the weather permits, 
which are the two more convenient locations sufficiently 
high above the sea surface. This is important, because 
any ship perturbations (shadow or reflections) are mini- 
mized when increasing the vertical distance between the 
operator and the sea surface, and then the horizontal dis- 
tance between the spot measured on the sea surface and 
the ship’s hull. 

The execution of a full sequence of measurements is as 
follows: 

1. Three dark current recordings (10s each); 

2. Three sun- viewing recordings (10 s each) for deriva- 
tion of the aerosol optical thickness; 

3. Three or four sea- viewing recordings (10s each); 

4. Three sun- viewing recordings (10 s each) for deriva- 
tion of the aerosol optical thickness; and 

5. Three dark current recordings (10s each). 

6.3 Coastal Sensors 

An automatic sun-tracking photometer model CE~318 
constructed by CIMEL Electronique is installed at the Cap 
Ferrat coastal site (Fig. 1). It has five standard Liters for 
440, 675, 870, 940, and 1,020 nnr and a polarized Liter for 
870 nnr. The instrument activates every 15nrin between 
0200 and 1400, but the exact measurement scheduling is 
preprogrammed by AERONET. The program consists of 
basic direct sun measurements — triplet observation, stan- 
dard measurement and Langley (plus basic sky measure- 
ments), Langley sky, alnrucantar, polarization, and prin- 
cipal plane. The procedures for data collection are as de- 
scribed on the AERONET Web sitef. 

6.4 Radiometer Characterization 

Characterization and calibration of held radiometers 
are crucial elements of a calibration and validation pro- 
gram for ocean color remote sensing satellites. The activi- 
ties performed in this segment are outlined in the following 
three subsections. 


f Full details are available at http : / /aeronet . gsf c . nasa . gov. 
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6.4.1 Absolute Calibration 


6.4.2. 1 The SQM- II 


The absolute calibration of the SPMR and SMSR with 
respect to NIST-traceable standards is performed every 
six months in the Satlantic optics calibration laboratory 
(Halifax, Canada). The objective of the absolute calibra- 
tions is to find a coefficient for each sensor to relate the 
raw instrument voltages (in digital counts) to the actual 
physical quantity of the parameter being measured in that 
environment. The method for doing this for radiance and 
irradiance sensors varies slightly and has been documented 
by Hooker et al. (2002a). 

For a radiance sensor (identified by Sid), the calibra- 
tion coefficient is computed using a plaque (identified by 
Tid) with a calibrated reflectance, along with a stan- 
dard lamp (identified by Aid) that has a calibrated irra- 
diance, Ej^ ( A). As a general procedure, the lamp is re- 
quired to be positioned on axis and normal to the center of 
the plaque at a distance d. Dark digital voltage levels are 
recorded with the radiance sensor capped and an average 
dark level, Asi D (A), is taken from these dark samples. 

The radiance sensor position allows a 45° view of the 
plaque with respect to the lamp illumination axis. Once 
the lamp is powered on, the voltage levels of each of the 
individual sensor channels are recorded. From these, an 
average calibration voltage for each channel, U 5 id (A), is 
obtained. The calibration coefficient is calculated as: 


^o(A) 


AR^(A)E^(A,50) 50 
H Sid (A) - D Sid (A) d 


( 1 ) 


The SQM-II is a portable, stable light source based on 
the design principles established by NASA and NIST with 
the original SQM (Johnson et al. 1998). The SQM-II was 
refined and built by Satlantic, Inc. The SQM-II consists 
of a lamp housing and a power source (deck box), which 
are connected by a 5m cable. A serial port provides the 
capability of monitoring and controlling the system with a 
PC. An internal memory, a liquid crystal display (LCD), 
and several buttons on the back allow manual control and 
monitoring. Several status light-emitting diodes (LEDs) 
are also used to provide a visual determination of the op- 
erational state of the instrument. 

The SQM-II has two sets of eight bulbs each. There is 
no individual bulb control, a set is either completely turned 
on or off. The light output of the second set, called the 
high power bank, is about three times brighter than the 
light output of the first set, named the low power bank. 
They are often referred to as HiBank and LoBank, respec- 
tively. These two lamp banks provide three illumination 
levels: low illumination when only LoBank is on, medium 
illumination when only HiBank is on, and high illumina- 
tion when both lamp banks are on. 

A sensitive internal detector at 490 nrn is used to mon- 
itor the stability of the lamps. The output from this de- 
tector and other internal sensors measuring temperatures, 
voltages and currents are monitored by a 20-bit analog-to- 
digital converter. The SQM-II setup is described in the 
SQM-II User’s Manual. 


where, d is given in centimeters. 

For an irradiance sensor (identified by Sid), the cali- 
bration coefficients are computed using an FEL standard 
lamp (Lid) with a calibrated irradiance 50). As 

a general procedure, the lamp is required to be on axis 
and normal to the face plate of the irradiance sensor at a 
distance, d. Similar to the radiance sensor, dark voltage 
levels (digital counts) are recorded with the radiance sen- 
sor capped and an average dark level, As id (A), is obtained 
from these dark samples. 

With the lamp powered on, the voltage levels of each of 
the individual sensor channels are recorded. From these, 
an average calibration voltage for each channel, Vs ID (A), is 
obtained. The calibration coefficient is calculated as: 


C* d (A) 


^(A,50) 
% d (A) - As id (A) 


r 1 2 

50 

~d 


( 2 ) 


where, d is given in centimeters. 


6.4.2 Stability 

The methodology and equipment used for tracking the 
stability of the main field radiometers is presented in the 
following three subsections. 


6. 4 .2. 2 Methodology 

To check the stability of the different radiometers in the 
field and to monitor the performance of the SQM-II itself, 
a calibration evaluation and radiometric testing (CERT) 
session and a data acquisition sequence (DAS) have been 
defined following the procedure described in Mueller and 
Fargion (2002). A CERT session is a sequence of DAS 
events, which are executed following a prescribed method- 
ology. Each DAS lasts for 3 min to represent enough data 
to statistically establish the characteristics of the instru- 
ment involved, referred to as a device under test (DUT). 
This records about 1,000 DUT records and 450 SQM-II 
records. A typical sequence for each CERT session is as 
follows: 

1. The equipment (lamp power supplies, the SQM-II 
fan and internal heater power supplies, lamp timers, 
etc.) is turned on at least 1 h before the CERT 
session begins. The total number of hours on each 
lamp set are tracked by recording the starting and 
the ending of hours on each lamp set. Each lamp 
bank functions separately, so the state of a given 
lamp bank is independent of the state of the other 
bank. 
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2. The SQM-II lamp bank will remain in this coarse 
stability state for 1 h, allowing enough time for the 
lamp and the system to thermally stabilize. It is 
not recommended that radiometer measurements 
be taken while in this state because the lamps have 
not yet reached their highest stability. After one 
hour, the lamp banks will automatically enter the 
fine stability state. 

3. If the mixture of radiometers used in the CERT 
session changes over time, at least one radiometer 
(preferably two of different types, i.e. , radiance and 
irradiance) should be used in all CERT sessions. 
This would be practically the case when the SPMR 
system and the buoy radiometers are not calibrated 
during the same CERT. 

4. S. Hooker advises to keep the banks in the fi n e 
stability state for another hour before starting the 
CERT, particularly in highly variable environments. 
This advice was followed even if the CERT took 
place in a laboratory. The warm-up period can 
be considered completed when the internal SQM-II 
monitor data are constant within 0.1%. The ra- 
diometric stability usually coincides with a thermal 
equilibrium as denoted by the internal thermistors. 

5. Inserting a DUT into the SQM-II produces a small 
reflective loading effect on the light chamber (i.e., 
some light is reflected back into the light chamber) , 
which affects what the radiometer detects. The 
reflectivity of a field radiometer, however, is con- 
stantly changing over time, because of the normal 
wear and tear associated with use, which alters the 
loading effect on the SQM-II monitor detector and 
the effective light field seen by the radiometer. 

6. Upon completion of the warm-up period, SQM-II 
monitor data are collected for the black glass (radi- 
ance) and white (irradiance) fiducials successively. 
A fiducial is a nonfunctional DUT, whose reflective 
surface must be carefully maintained over time so 
that its reflectivity remains essentially constant. It 
is also important that the position of the fiducial 
in the compartment is always the same and this is 
ensured by a fixed collar allowing only one position. 

7. Once these control measurements are completed, 
the individual radiometers (DUTs) are then tested 
sequentially. This begins with an OCR-200 that is 
designated specifically for use with the SQM-II and 
nothing else. Its use is restricted solely to SQM-II 
sessions to preserve the performance of its origi- 
nal calibration and, hence, serve as a verification 
of SQM-II stability between sessions. 

8. First, the DUT sensor head is inserted and secured 
in the SQM-II compartment. Data is collected over 
3 min. The DUT is secured using the three fasten- 
ers on the SQM-II, which clamp on a collar that 


has been fastened to the DUT. This collar is at a 
specific distance from the sensor head and of a par- 
ticular rotational position, determined by the flat 
side of the collar. This positioning of the flat part 
of the collar ensures that the DUT always enters 
the SQM-II in exactly the same way. 

9. The DUT is removed from the SQM-II slot and the 
caps are placed over the sensor heads to block out 
all possible light to the them. The black fiducial is 
then placed in the SQM-II and data is again col- 
lected from both SQM-II and DUT to provide a 
dark reading for the DUT and stability check for 
the SQM-II. Each time any file is recorded, the 
voltage at the SQM-II internal detector is noted. 

10. Once all instruments are tested, the SQM-II moni- 
tor data are collected a second time using the black, 
radiance, and irradiance fiducials successively. 

11. As soon as all measurements with a lamp bank are 
complete and its use is no longer required, the lamps 
are switched off to minimize aging and possible de- 
terioration. Where an alternative lamp bank com- 
bination is required for testing DUTs at a different 
light intensity, 2h warm-up time is once again re- 
quired to reach optimum stability. 

12. The practice of using all three light intensity levels 
was discontinued in 2003, because SQM-II perfor- 
mance changes caused the internal detector of the 
SQM-II to saturate during HiBank operation, so 
HiBank measurements were omitted. 

13. Before the SQM-II is finally shut down and the 
CERT session completed, after the lamps are pow- 
ered down, the ending number of hours on each 
lamp set is recorded. 

The SQM-II operation is precisely described step by step in 
the SQM-II User’s Manual edited by Satlantic, although 
some changes have been made to the method presented 
here to reduce the duration of each SQM-II session and to 
accommodate the aforementioned internal detector satu- 
ration during HiBank operation. 

6. 4- 2. 3 Measurement Schedule 

The objective of the scheduling of the SQM-II opera- 
tions is to monitor the variability and drift of the radiome- 
ters during the period of use between the (approximately) 
biannual absolute calibrations performed by Satlantic in 
Canada. The data collected in these SQM-II sessions aids 
in the reduction of uncertainties in the field data caused 
by changes in the performance of the sensors between cal- 
ibrations. When differences occur between calibrations, it 
cannot be assumed that this change has been a gradual 
and linear process during the period of activity. As soon 
as possible after the radiometers are received at LOV from 
Satlantic, they are tested on the SQM-II in the LOV optics 
laboratory using the method described in Sect. 6.2.2. 


22 



Antoine et al. 


Once the sensors have been used in the held, either 
fixed on the buoy or profiled from the ship, they are sub- 
jected to testing on the SQM-II. The exact scheduling of 
this testing depends on the duration of the held work. For 
the buoy radiometers, which are intended to be exchanged 
on a monthly basis, the SQM-II session is performed as 
soon as possible after they have been taken off the buoy 
and returned to the lab (once they have been wiped down 
with soapy water and a soft cloth). For the SPMR and 
SMSR, the SQM-II session is performed after each cruise. 

6.4.3 Immersion Factors 

The eighth SeaWiFS Intercalibration Round- Robin Ex- 
periment (SIRREX-8) conclusively demonstrated that im- 
mersion factors must be determined for each sensor if a 
less than 3% the uncertainty in irradiance determinations 
is to be achieved (Zibordi et al. 2002b). A water tank 
based on the Hooker and Zibordi (2005) design — the so 
called Compact Portable Advanced Characterization Tank 
(ComPACT) — was developed and is used for characteriz- 
ing buoy radiometers. 

The objective here is to determine specific sets of co- 
efficients for the radiometers, and, more importantly, to 
check whether these coefficients are changing over time 
because of their prolonged deployments at sea. For the 
moment, the results of these experiments do not show 
appreciable changes, although confirmation by additional 
measurements is still pending. More definitive conclusions 
also depend on future improvements in the reproducibility 
of the experiments, in order to meet the requested perfor- 
mance levels (Hooker and Zibordi 2005). 

6.5 Bio-Fouling 

Bio-fouling is a critical aspect of mooring applications. 
The development of new systems for preventing bio-fouling 
(chemicals or mechanisms) has not been considered in the 
project, at least up to now. Minimizing bio- fouling on 
instruments first requires that the structure they are in- 
stalled on is kept clean. When the mounting structure is 
bio-fouled, algal growth restarts more rapidly after the in- 
struments are cleaned, so an effective antifouling paint on 
the buoy superstructure is required. 

Concerning the radiometers, the quality of their mea- 
surements is preserved by the following ways: 

■ Copper sheets have been added around the radiome- 
ters’ bodies (since June 2005). 

■ The sensors are cleaned during each monthly main- 
tenance visit by divers equipped with brushes (since 
July 2006, the cleaning has been performed approx- 
imately every two weeks). 

■ Diffusers or filters are replaced if the absolute cal- 
ibration indicates it is timely to do so. This oper- 
ation can be performed after the upper superstruc- 
ture rotations, which occur about twice a year. 


■ Glass windows of the radiance sensors are kept as 
clean and smooth as possible, which is the best an- 
tifouling procedure following McLean et al. (1997). 

The above procedures are considered sufficient, at least 
as a first approach, especially given the oligotrophic-to- 
mesotrophic character of the mooring site. The cleaning 
of the radiometers every two weeks is the central element 
in the strategy to eliminate bio-fouling effects from the 
measurements of radiometric quantities. Concerning the 
transmissometers, copper rings have been installed around 
the emission and reception windows. A copper face plate 
has been installed onto the Hydroscat-II. 

All the above devices and procedures have proven to be 
efficient but cannot totally eliminate bio-fouling, so that 
corrections for possible residual effects will have to be ap- 
plied. The total effect of bio-fouling, which can be de- 
termined from the comparison of the measurements taken 
just before and just after cleaning of the sensors, will have 
to be distributed over time in some manner. 

7. DATA ANALYSIS 

Data analyses procedures are summarized in the follow- 
ing six subsections. This part is not meant to be exhaus- 
tive; only the analyses considered as the most important 
are presented. 

7.1 HPLC 

The methods described here were evaluated as part 
of three international SeaWiFS HPLC Analysis Round- 
Robin Experiment (SeaHARRE) exercises: SeaHARRE-1 
(Hooker et al. 2000), SeaHARRE-2 (Hooker et al. 2005) 
and SeaHARRE-3 (Hooker et al. 2007). The round robins 
were based on both natural samples and standard mixtures 
representing a wide range of Chi a concentrations from 
the hyper-oligotrophic conditions of the Mediterranean Sea 
and South Pacific Gyre to the very eutrophic waters of the 
Benguela, Moroccan, and Chilean upwelling systems. Four 
to eight laboratories from seven countries have taken part 
in these activities. 

As a result of these exercises, the LOV has demon- 
strated that satisfactory results are provided for the anal- 
ysis of Chi a and accessory pigments in seawater samples 
from different concentration regimes: the method used un- 
dergoes detailed quality control procedures, and a series of 
validation parameters are regularly monitored in order to 
ensure a good continuity in the data. 

7.1.1 Sample Collection 

Seawater samples were collected from Niskin bottles 
and filtered through 25 mm GF /F Whatman Liters (0.7 pm 
porosity). In most cases, 2.8 L were filtered for each sam- 
ple and the filters were placed in Petri slides, wrapped in 
aluminium and stored first in liquid nitrogen on board and 
then in a — 80° C freezer in the laboratory. 


23 



BOUSSOLE: A Joint CNRS-INSU, ESA, CNES, and NASA Ocean Color Calibration and Validation Activity 


7.1.2 Extraction Procedure 

Each filter is placed in 3 rnL of methanol (HPLC grade) 
containing an internal standard! • After 30 min at — 20°C, 
the filters are disrupted by ultrasonication using an ultra- 
sonic probe and returned to the freezer. Another 30 min 
later, the sample is clarified through a 25 mm GF/C What- 
man filter (1.2 pm porosity). The filtrate is finally stored 
at — 20°C until analysis (within 24 h). 

7.1.3 HPLC Analysis 

All parts of the HPLC system at the LOV are Agilent 
Technologies products: 

• A degasser (1100 model); 

• A binary pump (1100 model); 

• An autosampler (1100 model) with built-in temper- 
ature control (4°C) and automatic injection (Rheo- 
dyne valve) for mixing the sample with the ammo- 
nium acetate (1 N) buffer; 

• A diode array detector (1100 model) with measure- 
ments at 222 nrn (for the vitamin E acetate inter- 
nal standard), 440 nrn (for carotenoids), and 667 nrn 
(for chlorophylls and degradation products); and 

• A fluorescence detector with excitation and emis- 
sion wavelengths at 417 and 670 nm, respectively. 

Two successive analytical methods have been used for 
pigment analysis. The first method is based on a gradient 
separation between a methanobammonium acetate (70:30) 
mixture and a 100% methanol solution, comparable to that 
described by Vidussi et al. (2000), but with a few differ- 
ences allowing for improvement of sensitivity and peak res- 
olution. Modifications to this method to separate certain 
peaks and increase sensitivity included a) a flow rate of 

0.5 mL min , and b) a reversed phase chromatographic 
Cg column with a 3 mm internal diameter (Hypersil MOS 
3 pm). The gradient used is presented in Table 3. Col- 
umn temperature is maintained at 25° C and the injection 
volume is 200 pL. 

The second method is also a gradient separation be- 
tween a methanoktetrabutyl-amonium acetate (70:30) mix- 
ture and a 100% methanol solution, comparable to that de- 
scribed by Van Heukelem and Thomas (2001), except for a 
few differences allowing improvement of the sensitivity and 
peak resolution. Modifications to this method include a) a 
flow rate of 0.55 mL min -1 , b) a reversed phase chromato- 
graphic Cg column with a 3 mm internal diameter (Zorbax 
Eclipse XDB, 3.5 pm). The gradient used is presented in 
Table 3. Column temperature is maintained at 60°C and 
the injection volume is 250 pL. 


f Samples before March 2003 were analyzed with trans-(3- apo- 
B'-carotenal as the internal standard, which was replaced 
with vitamin E acetate. 


Table 3. The gradient used with the HPLC col- 
umn. The time is in minutes, and the percentages of 
solvents A and B are given in the last two columns. 


Step 

Time 

A [%] 

B [%] 

1, Start Method 1 

0 

80 

20 

2 

4 

50 

50 

3 

18 

0 

100 

4, End Method 1 

22 

0 

100 

1, Start Method 2 

0 

90 

10 

2 

25 

5 

95 

3, End Method 2 

28 

5 

95 


The different pigments that are quantitated are pre- 
sented in Table 4. 

7.1.4 Calibration 

Response factors for eight pigment standards provided 
by DHI Water and Environment (Hprsholm, Denmark) are 
determined by spectrophotometry (Perkin Elmer) followed 
by HPLC analysis: 

1. Peridinin, 

2. 1%-Butanoyloxyfucoxanthin, 

3. Fucoxanthin, 

4. lSy-Hexanoyloxyfucoxanthin, 

5. Alloxanthin, 

6. Zeaxanthin, 

7. Chlorophyll b, and 

8. Chlorophyll a. 

These response factors are then derived to compute the 
specific absorption coefficients of these eight major pig- 
ments for the HPLC system. 

Specific absorption coefficients for divinyl chlorophyll a 
and divinyl chlorophyll b are computed knowing a) the 
specific extinction coefficients of chlorophyll a and chloro- 
phyll 6; b) the measurement of the absorption of chloro- 
phyll a and divinyl chlorophyll a (or chlorophyll b and di- 
vinyl chlorophyll b) at 440 nm when the spectra of both 
pigments are normalized at their red maxima; and c) that 
both pigments are considered to have the same molar ab- 
sorption coefficient at this red maxima. For the remaining 
pigments, their specific extinction coefficients were either 
derived from previous calibrations or from the literature 
(Jeffrey and Vesk 1997). 

7.1.5 Quantitation 

The Agilent Technologies Chemstation software is used 
for conducting the analysis, as well as for post-analysis pro- 
cessing. This includes peak integration and spectral iden- 
tification. Peak identification is manually verified by re- 
tention time comparison and observation of the absorption 
spectra. Quantification is based on peak area related to the 
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Table 4. The list of pigments detected by HPLC at the LOV, their detection wavelengths, and their possible 
coelution with another pigment. The variable forms, which are used to indicate the concentration of the 
pigment or pigment association, are patterned after the nomenclature established by the Scientific Committee 
for Oceanographic Research (SCOR) Working Group 78 (Jeffrey and Vesk 1997). Abbreviated forms for the 
pigments are shown in parentheses. 


Variable 

Pigment 

Detection 

Coelution 

[Chi c 3 ] 

Chlorophyll C 3 (Chi C 3 ) 

440 [nm] 


[Chlide a] 

Chlorophyllide a (Chlide a) 

667 

Coelutes with Chlorophyll ci+c 2 

[Chi c 1+ c 2 ] 

Chlorophyll ci + c 2 ( Chi ci + c 2 ) 

440 


[Phide a] 

Phaeophorbide a (Phide a) 

667 

Coelutes with Peridinin 

[Peri] 

Peridinin (Perid) 

440 


[But] 

1 9 B ut anoy loxyfucoxant hin ( B ut- fuco ) 

440 


[Fuco] 

Fucoxanthin (Fuco) 

440 


[Hex] 

1 9'-Hexanoyloxyfucoxanthin (Hex-fuco) 

440 


[Neo+Vio] 

Neoxanthin and Violaxanthin (Neo+Viola) 

440 

Coelution Pair 

[Diad] 

Diadinoxanthin (Diadino) 

440 


[Alio] 

Alloxanthin (Alio) 

440 


[Diato] 

Diatoxanthin (Diato) 

440 


[Zea] 

Zeaxanthin (Zea) 

440 


[Lut] 

Lutein (Lut) 

440 


[Chi ci] 

Nonpolar chlorophyll ci (Chlci) 

440 


[TChl b] 

Total chlorophyll b f (TChl b) 

440, 667 

Coelution Pair 

[Croco] 

Crocoxanthin (Croco) 

440 


[DVChl a] 

Divinyl chlorophyll a (DVChl a) 

440 


[Chi a] 

Chlorophyll a\ (Chi a) 

440, 667 


[Chi c 2 ] 

Nonpolar chlorophyll c 2 (Chi c 2 ) 

440 


[Caro] 

Carotenes§ (Caro) 

440 

Coelution Pair 

[Phytin a] 

Phaeophytin a (Phytin a) 

667 



f The sum of divinyl chlorophyll b and monovinyl chlorophyll b. 
i The sum of monovinyl chlorophyll a plus allomers and epimers. 
§ The sum of /3/3- Car plus /3e-Car. 


specific extinction coefficient and concentrations are given 
in milligrams per cubic meter. When two pigments tend 
to coelute, their identification is first done spectrally then 
they are summed, e.g., chlorophyll C 1 +C 2 , or total chloro- 
phyll b (chlorophyll b plus divinyl chlorophyll b ) . 

7.1.6 Validation 

A solution of methanol containing the internal standard 
is injected every 12 samples during the analytical sequence. 
The average of these control analyses provides the refer- 
ence peak area for the internal standard. The standard 
deviation of these analyses provides information about the 
precision, as well as the stability of the instrument dur- 
ing the analytical sequence. Total Chi a (TChl a) values 
are compared to the fluorescence signal from CTD bot- 
tle data in order to detect possible inaccuracies. TChl a 


values are also compared to particulate absorption mea- 
surements, which have been carried out on the same filters 
just before pigment extraction (Sect. 7.2). 

7.1.7 Method Performance 

This HPLC analytical method has proved to be par- 
ticularly sensitive, with detection limits of approximately 
0.0004 mg m -3 for [Chi a] and a resolution between divinyl 
chlorophyll a and chlorophyll a greater than 1.5, which al- 
low it to be particularly adapted to the oligotrophic waters 
of the Mediterranean. The precision of the method is char- 
acterized by a variation coefficient of about 0.4%. 

7.2 Filter Pad Absorption 

The following two subsections describe the measure- 
ment of absorption by particles, which is performed on the 
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same filter that is subsequently used for the HPLC pigment 
analysis previously described. 

7.2.1 Laboratory Analyses 

The protocol for analyzing the spectral absorption co- 
efficients of particulates is based on the NASA Ocean Op- 
tics Protocols for Satellite Ocean Color Sensor Validation 
(Mitchell et al. 2003). Sampling is performed in the same 
way as for the HPLC pigments as described in Sect. 8.1 
where the seawater samples are collected from Niskin bot- 
tles and filtered through 25 mm GF /F Whatman Liters 
(0.7 pm porosity). In most cases, 2.8 L are filtered for each 
sample and the Liters placed in petri slides, wrapped in 
aluminum foil and stored Lrst in liquid nitrogen on board, 
and then in a — 80° C freezer in the laboratory. The sample 
Liters used for the A p analysis are also used afterwards for 
HPLC analysis. This method has been tested to ensure 
that using the same Liter for the two methods does not 
compromise either measure. 

The spectrophotometer used is a Perkin Elmer Lambda 
19 (L19) dual-path spectrometer with an integrating sphere 
compartment attached. The instrument scan range is set 
from 750-350 nrn with a scan interval of 1 nrn. After a 1 h 
warm-up period, the machine is ready to use. The samples 
are taken out of the — 80°C freezer and stored temporarily 
in an ice-cooled cold box. Each sample, in its petri dish, is 
then placed in a dark box on the bench to defrost for 5 min. 
This defrosting period is kept to a minimum because the 
pigment composition once thawed becomes susceptible to 
degradation which could effect the analyses, particularly 
the HPLC. 

The L19, being a dual-path spectrometer, has a slot 
for a reference material and another for the sample. These 
slots are actually the outer walls of the integrating sphere 
at the two positions where the two beams enter the sphere. 
For the reference path, rather than use a blank Liter, which 
can create variations in the spectra attributable to differ- 
ences in the level of nroistness of the Liter, the pathway 
is left open (without any Liter blank), thus using an air 
blank. In the sample material slot, a blank GF/F Liter 
is used. This blank has been soaked in distilled water for 
12 h and had excess water drained off before being placed 
in the sample slot. 

An autozero is performed with the blank Liter in the 
sample path and the reference path open, which sets the 
absorbance values measured at each wavelength to zero. 
As a veriLcation of this and to provide a baseline, the blank 
is then scanned to measure absorbance. The spectra pro- 
duced from this should be Lat and close to zero (within 
0.005). If this is not the case, the autozero is repeated and 
the blanks scanned again. 

Once satisfactory baselines are achieved, then the sam- 
ple Liters can be scanned. For both cases, two scans are 
made, which should produce two similar replicates. If sim- 
ilarity is not found, the Liter position is checked and the 


sample scanned again. Drying of the Liter between repli- 
cates may cause some vertical offsetting, which is accept- 
able if the spectral shape for each replicate is similar. 

The spectra are saved as a Lie in American Standard 
Code for Information Interchange (ASCII) format with two 
columns: wavelength and absorbance. The mean of the ab- 
sorbance values of replicates for each sample is calculated 
and any offset from zero removed by subtracting the ab- 
sorbance value at 750 nrn from the entire spectrum. 

7.2.2 Data Postprocessing 

The optical densities are then corrected for any nonzero 
signal in the NIR (750 nrn), and transformed into absorp- 
tion coefficients by accounting for the so-called f3 effect. 
The total particulate absorption spectra are then numeri- 
cally decomposed into a phytoplankton and a detritus ab- 
sorption spectra following Bricaud and Stramski (1990). 

7.3 Total Suspended Matter (TSM) 

TSM is determined with a gravimetric method based on 
the protocol described by Strickland and Parsons (1972). 
Before sampling, 25 mm diameter Whatman GF /F glass 
Lber Liters (0.7 pm particle retention size) are pre-rinsed 
in methanol, pre-combusted at 450° C for 4h, and pre- 
weighed to six signiLcant digits with a Metier- Toledo elec- 
trobalance (ambient temperature and humidity are also 
recorded). 

In the held, 1.5-5.6L seawater samples (the chosen vol- 
ume depends on the amount of material in the water) from 
below the surface are collected in triplicate using vacuum 
Lltration. The Liters and Llter-rims are rinsed with 20- 
30mL of milli-Q water in order to eliminate the presence 
of salt crystals. They are placed in polycarbonate petri 
dishes and stored in liquid nitrogen and at — 80°C until lab- 
oratory analysis. The Liters are dried in an oven at 75°C 
for 24 h and placed in a dessicator before being weighed to 
six digits of precision (ambient temperature and humid- 
ity are recorded). The TSM is calculated as the difference 
between the sample weight minus the weight before Lltra- 
tion (in grams), and then divided by the Lltration volume 
in cubic meters (concentrations are, therefore, grams per 
cubic meter). 

7.4 AC-9 and ECOBB3 Data Processing 

When a CTD cast is performed, the CTD and AC-9+ 
start times are not synchronized. In addition, each in- 
strument has a speciLc time constant and the water cir- 
culation in the tubing generates additional temporal off- 
sets. In order to account for these time lags, the verti- 
cal proLles for absorption and attenuation must be shifted 
in order to match the vertical proLles of other data such 
as the chlorophyll Luorescence and the temperature. The 
matching is either performed manually using clearly iden- 
tihable features in the proLles (i.e. , peaks) or automati- 
cally when the time lag between the start of the various 
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instruments is accurately known. Once the vertical shift 
has been determined for the AC-9+ data, it is also applied 
to the ECO-BB3 data. 

When all data have been matched in the vertical, the 
data processing for the AC-9+ essentially follows the pro- 
tocol described in the WETLabs AC- 9+ User’s Guide , in 
terms of temperature, salinity, and scattering corrections 
(for the latter, the so-called method 3 is used). For the 
ECO-BB3, the data processing follows the protocols de- 
scribed in the Wetlabs Scattering Meter, ECO-BB User’s 
Guide. 

For the AC-9+ data, temporal drifts in the absolute 
calibration of the device (from cruise to cruise) are cor- 
rected as described in Oubelkheir et al. (2005), i.e. , av- 
erage absorption and attenuation spectra are determined 
over depths from 350-400 m immediately after an abso- 
lute calibration. These spectra are then used as reference 
spectra, i.e., all spectra subsequently measured at 400m 
are constrained to be equal to the reference spectra. Any 
subsequent deviation of the coefficients measured at 400 m 
from these references is attributed to an instrument drift. 
The deviation, if any, is simply removed from the entire 
profile. 

At the end of these processing steps, a verification of 
the AC-9+ absorption data is performed by comparing 
the values at 676 nrn to the phytoplankton absorption at 
the same wavelength, as derived from discrete samples 
(Sect. 7.2). Attenuation at 660 nrn is also compared to 
the beam attenuation coefficient determined by the C-Star 
transmissometer (installed on the CTD package since the 
beginning of 2006). 

7.5 AOPs from SPMR Profiles 

The SPMR is a central element in the BOUSSOLE ac- 
tivity. It provides reference data against which the data 
collected from the buoy can be checked, and it also pro- 
vides water column information that is unobtainable from 
the buoy (near-surface) measurements. 


where z is depth, E u (z 0 ) is the shallowest value of E u (z) for 
which the tilt is less than 2° , and K u is the diffuse attenu- 
ation coefficient for the upward irradiance computed from 
the measurements of E u (z) collected at different depths 
between z 0 and z 0 + 20 m. 

Several interpolation procedures designed to derive the 
E u ( 0 _ ) from the vertical profile of E u (z) have been tested 
against true values of E u ( 0”) (i.e., values directly mea- 
sured below the sea surface by installing the radiometer on 
a floating frame), and the method that eventually provided 
the closest values to the true ones was selected. This ex- 
perimental work, which is not further detailed here, is just 
mentioned to indicate that the contribution to the overall 
uncertainty budget of the extrapolation error is tentatively 
minimized. 

The above- water solar reference data, E s , are corrected 
to account for the loss of irradiance at the air-sea interface: 

* W 0 -) = e-y Gjir < 4 > 

where the mean transmission of the sea surface for sky and 
Sun irradiance, expressed by 1 — p, is equal to 0.957 (±3% 
according to atmospheric turbidity and solar elevation), 
and the internal reflectance (accounted for by 1—fR, where 
is f = 0.489), varies slightly with R. With a mean R value 
of 3% this term is equal to 0.985 (±1.5% if R varies between 
0-6%). Assuming these two terms are constant, 

E d ( 0-) = 0.97 E s , (5) 

where the value of E s is obtained from the first 10 s of 
recording starting after the release of the SPMR (this cor- 
responds approximately to the upper 5m of the descent), 
to which a fit is adjusted in order to eliminate variations 
in E a , which are only due to the tilt of the SMSR (which 
is not installed on a girnbal). This technique provides sim- 
ilar results as compared to just picking the measurements 
taken for tilt angles less than 1°. 


7.5.1 Measurement Suite 

What is measured in the water column is the upward 
irradiance, E u (z, A), at 13 wavelengths from 412-865 nrn 
(Table 1), plus the downward irradiance at the same wave- 
lengths, E d (z, A). The latter is not used in the computa- 
tion of the reflectance. The above-water downward irradi- 
ance, Eg, which is often referred to as a deck cell measure- 
ment, because is recorded on deck (in this case at the bow 
of the ship) and again at the same 13 wavelengths. 

7.5.2 Corrections and Extrapolations 

From the vertical profile of E u (z, A), the upward irra- 
diance at null depth (0"), i.e., immediately below the sea 
surface, is obtained as (A is omitted for brevity): 

E u ( 0-) = E u {z 0 )e zK \ (3) 


7.5.3 Computed Parameters 

The reflectance R is then 


E,(0-) 

E d { 0-)’ 


( 6 ) 


remembering that before the above ratio is formed, E u { 0 _ ) 
is corrected for instrument self shading as per Gordon and 
Ding (1992). In this correction, the instrument radius, 
the total absorption coefficient, and the ratio between the 
direct-sun and diffuse-sky irradiances (r^) are taken into 
account. All of the Satlantic sensors used here have a com- 
mon radius (4.5 cm), the absorption is computed following 
Morel and Maritorena (2001) using the measured Chi a 
concentration, and r d is computed following Gregg and 
Carder (1990). 
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An average diffuse attenuation coefficient for the down- 
ward irradiance in the upper layers is computed as 

log [E d {z)/E d { 0-)] 

K d = , {<) 

where z is taken equal to be 20 m to ensure that the per- 
turbations of the light field caused by surface waves are 
not introducing errors into the computation. 

7.6 Buoy Data Preprocessing 

The binary log hies for the instruments connected to 
the DATA-100s are converted by a program called Satcon 
into ASCII format hies. The Satcon program uses the 
calibration coefficients and format information in the cali- 
bration (* . cal) hies, created by Satlantic during the most 
recent absolute calibration, to convert the raw counts to 
physical units. The compass and tilt sensor data are con- 
verted by Satcon using coefficients and information in a 
text definition file (*.tdf), provided by Satlantic. The 
ASCII hies produced by Satcon are given the hie exten- 
sion * . dat . 

The Hydroscat-II log hie data are converted to physical 
units using HOBI Labs’ Hydrosoft software and the HOBI 
Labs calibration hies and output to (*.dat ASCII) hies. 
There are no valid time values associated with these hies, 
except the time of hie creation in the filename. 

The CTD * . log hies do not require preprocessing be- 
fore they are concatenated and merged with the other in- 
strument data. There is no time data associated with these 
hies except for the hie creation time in the log hie name. 

A Matlab processing script removes header lines from 
all the ASCII * . dat hies then, where multiple hies for one 
day exist, these are concatenated into a single hie for each 
day and for each instrument. In order to have a single 
day hie for all the instruments, the timeframe of the 4 m 
DATA-100 hie is used as a standard for each of the other 
hies. The number of rows in the hnal hie, therefore, is 
always the same as the 4 m DATA-100 hie. Two processes 
are required for integrating the other instrument hies, one 
for the hies that have their own time stamps already and 
one for those that do not. The latter relates to the CTD 
and Hydroscat-II hies. 

Where a time stamp is present for an instrument, inte- 
grating the data time frame for this instrument into the 4 nr 
DATA-100 timeframe is performed by hnding the closest 
time stamp in the 4 m DATA-100 hie that is less than the 
time stamp for each line of the hie of interest and within 
the same 1 min sampling period. 

For the Hydroscat-II and CTD data, which are not 
time stamped, their integration into the 4 nr DATA-100 
time base is performed using the proportion of each 1 nrin 
4 nr DATA-100 sampling period to the entire day of 4 nr 
DATA-100 data. The progressive proportions of each of 
these sampling periods for a day of data is then used to 


divide up the day of CTD and Hydroscat-II data to the 
same ratio. Once these hies are divided into theoretical 
1 nrin slots, the data lines within each slot are then dis- 
tributed evenly across the true 1 nrin sampling period of 
the 4 nr DATA-100. 

For the seven radiometers (OCR-200 and OCT200), a 
dark correction is performed by subtracting, from each ra- 
diometer data series, the mean data values collected be- 
tween the hours of midnight and 0200 hours for each in- 
strument. After processing, general radiometer health can 
be assessed for each day by monitoring these dark values. 
The dark values for each instrument are, therefore, saved 
to a hie where they can be easily checked for changes. The 
buoy data after preprocessing is presented as large (approx- 
imately 13 MB) individual hies containing all instrument 
data for a single day from midday to midnight. 

7.7 IOPs from the Buoy Data 

The processing of the huorescence and IOP data col- 
lected on the buoy simply follows the procedures recom- 
mended by the instrument manufacturers, i.e., the Chelsea 
User Guide for MINItracka II C In-Situ Fluorimeter HB 
176 for phytoplankton huorescence, the Wetlabs C-Star 
Transmissometer User Guide for the beam attenuation co- 
efficient at 660 nm, and the HOBILabs Hydroscat-II Spec- 
tral Backscattering Sensor User's Manual (Revision F, Oc- 
tober 25, 2004) for the backscattering coefficient. Concern- 
ing the latter, the recommendations of Mahone and Dana 
(1997) and Boss and Pegau (2001) are used to transform 
the measurement of the volume scattering function at 140° 
into the particle backscattering coefficient. 

7.8 AOPs from Buoy Data 

The apparent optical properties that are derived from 
the radiometry measurements collected by the buoy are 
providing the core of the satellite matchup data set. 

7.8.1 Measurement Suite 

The in-water optical measurements are made at two 
depths (nominally 4 and 9 m) and at seven wavelengths 
from 412-681 nm (Table 1). Three types of measurements 
are acquired: a) the upward irradiance, E u (z, A); b) the 
downward irradiance, EJ(z, A), which is not used in the 
computation of the reflectance; and c) the upwelled (nadir) 
radiance, L u (z,A). The above-water downward irradiance, 
E s ( A), is recorded on top of the buoy at approximately 
4.5m above the surface (at the same seven wavelengths). 

7.8.2 Computed Parameters 

The first step in the data processing scheme consists of 
a data reduction, which aims at producing one represen- 
tative value for each of the 1 min sequences of acquisition. 
During these sequences, about 360 measurements are col- 
lected (acquisition rate of 6 Hz), which may include some 
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outliers and are usually noisy because of the fluctuations 
of the radiative field due to the wave-roughened air-sea 
interface. Taking the median of the 360 measurements has 
been found to be a simple and robust way of filtering the 
data. 

From the two values of E u (z , A), the upward irradiance 
at null depth z = 0“ (immediately below the sea surface) 
is obtained as (omitting wavelength for brevity): 

E u ( 0-) = E u {A)e ZK \ (8) 

where Z is the measurement depth (not exactly 4 m when 
the buoy is lowered or when swell is going through the su- 
perstructure) , K u is the diffuse attenuation coefficient for 
the upward irradiance computed from the measurements 
of E u collected at the two depths. The same extrapolation 
can be performed from the measurement at 9 m. 

From the two values of L u ( A) at 4 and 9 m, the up- 
welling radiance at null depth is obtained as: 

L u (0~) = L U (A) e ZKh , (9) 

where Z is the measurement depth, and Kl is the dif- 
fuse attenuation coefficient for the upwelling radiance com- 
puted from the measurements of L u collected at the two 
depths. 

Before combining these values to the above- water refer- 
ence measurement, E s , in order to compute the reflectance, 
E s is corrected for the buoy tilt. The correction is a func- 
tion of the orientation of the two axes of the tilt measure- 
ment with respect to the solar azimuth angle, and com- 
putes the ratio of the diffuse (not affected by the tilt) to 
direct (affected) light for clear-sky conditions. 

The computation of R. follows (6) and Kd follows (7), 
where the depth z can be from either measurement depths 
(nominally 4 and 9 m). This coefficient is also computed 
between these two measurement depths. 

7.9 Sun Photometer Data 

The first data set that comes from the processing of the 
sun photometer data includes the aerosol optical thickness 
at five wavelengths (Table 1). Numerous methods have 
been developed in the past years, which are dedicated to 
the retrieval of information about aerosol types, usually 
particle size distributions, refractive indices, or phase func- 
tions from measurements of sky radiances collected either 
in the principal plane, following almucantars, or in the so- 
lar aureole. Good examples of these activities are provided 
by Santer and Martiny (2001), Dubovik and King (2000), 
Duboviket al. (2000), and Nakajimaet al. (1983 and 1996). 

A similar method has been developed, which is distin- 
guished by using the degree of polarization, as measured 
at 870 nrn by the sky photometer. Using this additional 
piece of information, in principle, decreases the ambigui- 
ties of the inversion process. In order to find the best can- 
didate aerosol model that allows a reconstruction of the 


sky radiances distribution observed in the principal plane, 
the method either uses the trial-and-error principle or a 
more sophisticated inversion algorithm based on the use 
of a neural network. A radiative transfer code is used in 
these inversions, which is based on the successive orders of 
scattering method and uses the vector theory Ordres Suc- 
cessifs Ocean et Atmosphere , (OSOA) code (Deuze et al. 
1989 and Charni et al. 2001). 

In regards to the neural network approach, which has 
the advantage of allowing fast processing of time series 
and has the ability of properly dealing with measurement 
errors, the important step is the training of the neural net- 
work. This training has also been performed by using the 
OSOA radiative transfer code, and using aerosol models 
with varying indices of refraction and particle size distri- 
butions following either a Jiinge model or a log-normal 
distribution with two modes (small and large particles). 
A realistic noise that accounts for instrumental errors is 
added to the synthetic data. 

8. INITIAL RESULTS 

In the following sections, sample results are provided 
for the three basic data collection opportunities: the ship 
campaigns, the buoy daily acquisitions, and the sun pho- 
tometer data. The purposes of the presentations are to 
establish the extent of the sampling program, in terms of 
the type and quality of data products being produced, and 
to illustrate the site characteristics. 

8.1 Data Collection Summary 

The following three sections present a quantitative sum- 
mary of the data collected from the ship, from the buoy, 
and at the coastal AERONET site. 

8.1.1 Ship Data 

Table 5 provides a list of monthly maintenance cruises 
until April 2006 (and continuing since then on a monthly 
basis), and indicates how many data collection events have 
been executed for each instrument category. Data collec- 
tion is either in the form of a vertical profile of an appropri- 
ate portion of the water column (SPMR, CTD, and IOPs), 
individual temporal sequences based on a proper amount 
of time or number of samples (CIMEL and SIMBADA), or 
discrete samples of filtered seawater taken from a variety 
of depth levels {A p , HPLC, and TSM). 

8.1.2 Buoy Data 

Table 6 indicates the number of buoy deployment days 
for the successive deployment sequences. The gap in 2001 
is associated with the redesign effort of the original (ver- 
sion 1) buoy to correct some deficiencies noted in the first 
deployment (but note the continuation of the site charac- 
terization during this time period as shown in Table 5). 
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Table 5. A list of the monthly maintenance cruises from 2001-2005 (inclusive) and the types of measurements collected. 
The cruise identification (ID) number provides a sequential enumeration of the data collection sequence. All scheduled 
cruises are shown to establish the level of effort undertaken. If a cruise was cancelled (e.g., ID 27), perhaps because 
of bad weather or a scheduling conflict, the number of days is shown as zero and there are no data entires for this 
maintenance cycle. A checkmark indicates data were collected and a numerical value indicates how many discrete 
samples, profiles, or acquisition sequences were acquired. 


Year 


Dates 


: ID 

Days 

SPMR 

CTD 

IOPs 

CIMEL 

SIMBADA 

Ap 

HPL( 

24 Jul 

1 

3 

19 




8 


9 

4 Sep 

2 

3 

21 

6 



20 



4 Oct 

3 

2 

2 






20 

5 Nov 

4 

1 

14 

3 





2 

4 Dec 

5 

2 

11 




7 


9 

14 Dec 

6 

0 








23 Jan 

7 

3 

6 

5 





13 

27 Feb 

8 

0 








24 Mar 

9 

2 

11 

2 



6 


9 

12 Apr 

10 

0 








22 May 

11 

3 

12 

2 





28 

28 Jun 

12 

4 

24 

10 



11 


31 

22 Jul 

13 

3 

26 

5 

/ 


10 


31 

4 Sep 

14 

3 

9 

3 

/ 


7 


20 

7 Oct 

15 

2 

11 

2 

/ 


3 


23 

28 Oct 

16 

0 








29 Nov 

17 

2 

14 

4 

/ 




29 

9 Dec 

18 

0 








9 Feb 

19 

2 

18 

4 

/ 


5 

21 

47 

19 Mar 

20 

2 

8 

1 

/ 


4 

20 

27 

7 Apr 

21 

1 

11 

2 

/ 


12 

20 

20 

29 May 

22 

3 

36 

13 

/ 


30 

14 

50 

29 Jun 

23 

3 

24 

10 

/ 


25 

28 

47 

15 Jul 

24 

3 

47 

14 

/ 


26 

57 

87 

31 Aug 

25 

0 








25 Sep 

26 

3 

10 

9 

/ 


1 

29 

36 

19 Oct 

27 

0 


1 






7 Dec 

28 

3 

12 

10 

/ 



32 

48 

18 Dec 

29 

2 

12 

7 




19 

37 

1 Feb 

30 

5 

20 

13 

/ 


5 

32 

39 

7 Mar 

31 

3 

15 

9 

/ 


5 

20 

32 

12 Apr 

32 

4 

11 

10 

/ 



27 

36 

4 Jun 

33 

4 

30 

12 

/ 



28 

34 

7 Jul 

34 

3 

11 

10 

/ 



10 

43 

17 Aug 

35f 

9 

54 

92 

/ 

105 

180 

247 

279 

23 Oct 

36 

3 

5 

1 

/ 

6 


4 

14 

24 Nov 

37 

3 

9 

11 

/ 



12 

49 

4 Feb 

38 

3 

13 

13 

/ 



11 

43 

26 Feb 

39 

3 

5 

10 

/ 

8 


14 

28 

28 Mar 

40 

3 

6 

11 

/ 

5 


26 

38 

18 Apr 

41 

4 

3 

8 

/ 

9 


27 

36 

28 May 

42 

3 

8 

8 

/ 

12 


21 

21 

29 Jun 

43 

3 

9 

12 

/ 

12 


25 

33 


2001 


2002 


2003 


2004 


2005 


22 Jul 
2 Sep 

2 Oct 

3 Nov 
2 Dec 

12 Dec 

21 Jan 

25 Feb 

22 Mar 
10 Apr 

20 May 

26 Jun 
20 Jul 

2 Sep 
5 Oct 

26 Oct 

27 Nov 
7 Dec 
7 Feb 

17 Mar 
5 Apr 
27 May 

27 Jun 

13 Jul 
29 Aug 

23 Sep 
16 Oct 

5 Dec 
16 Dec 

28 Jan 

3 Mar 
9 Apr 
2 Jun 
5 Jul 

31 Jul 
15 Oct 

21 Nov 
31 Jan 
23 Feb 

25 Mar 
15 Apr 

26 May 

27 Jun 


f Cruise 35 was the Advanced Optical Properties Experiment (AOPEX). It was carried out on a larger ship (R/V Le Suroit), 
and involved 17 scientists from LOV plus other European and American laboratories. Several additional parameters, such as 
the volume scattering function of particles, were collected, which are not indicated in this table. 
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Table 5 (cont.). 


A list of the monthly maintenance cruises for 2006 and the types of measurements collected. 


Year Dates 


2005 


2006 


27 Jul -»• 29 Jul 

9 Sep — > 11 Sep 
8 Oct — > 11 Oct 
5 Nov — > 8 Nov 

13 Dec — > 16 Dec 

1 Feb — > 4 Feb 

28 Feb — > 4 Mar 

14 Mar — > 17 Mar 

3 Apr — > 6 Apr 


ID 

Days 

SPMR 

CTD 

IOPs 

CIMEL SIMBADA 

Ap 

HPLC 

TSM 

44 

3 

7 

11 

/ 

5 

38 

51 


45 

3 


7 

/ 

2 

19 

21 


46 

4 

7 

9 

/ 

13 

22 

28 


47 

4 

15 

12 

/ 

4 

30 

53 

9 

48 

4 

2 

6 

/ 

2 

9 

10 

3 

49 

4 

6 

9 

/ 


23 

28 

4 

50 

1 


7 

/ 


11 

11 


51 

4 

14 

10 

/ 


20 

29 

6 

52 

3 

14 

10 

/ 

9 

31 

36 

9 


The main purpose of the redesign was to improve the per- 
centage of time for which the stability requirements — in 
terms of inclination — were satisfied. 


Table 6. A summary of the buoy data collected at 
the BOUSSOLE site. 


Year 


Dates 

ID 

Days 

2000 

20 Jul 

-► 20 Oct 

1 

90 

2002 

16 May 

— > 6 Jun 

2 

11 

2003 

6 Sep. 

— > 6 Dec 

3 

90 

2004 

4 Mar. 

— > 7 Apr 

4 

35 


11 May 

-»■ 23 Jul 

5 

72 


5 Aug. 

— ► 31 Dec 

6 

149 

2005 

18 Feb 

— > 2 Jun 

7 

105 


17 Jun. 

— > 31 Dec 

8 

197 

2006 

1 Jan 

— > 10 Jun 

8 

161 


11 Jun. 

— > To Date 

9 



8.1.3 Sun Photometer Data 

Table 7 presents how many days of data collection have 
been successfully completed at the sun photometer site. 


Table 7. A summary of the sun photometer data 
collected at the Cap Ferrat coastal site. 


Year 


Dates 

ID 

Days 

2002 

3 Jul 

— > 31 Dec 

1 

181 

2003 

1 Jan 

— > 15 Apr 

It 

105 

2004 

6 Jan 

— > 10 Nov 

2 

309 

2005 

16 Feb 

— > 5 Oct 

3 

231 


10 Feb 

— > 29 Mar 

4 

47 

2006 

26 May 

— ► To Date 

5 



f The completion of the first deployment sequence. 


8.2 Physical Conditions 

Figure 13 is a time series of the potential temperature 
over the upper 200 m of the water column spanning the 
time period from July 2001 to January 2006. The well- 
known winter mixing and spring stratification clearly ap- 


pear, as well as an increase of the summer surface tempera- 
ture from year to year. The summer of 2003 was one of the 
warmest recorded over the last several decades in northern 
Europe, and the signature of this unusual climate is seen 
here in the temperature record. 

Another noteable feature in Fig. 13 is seen in the win- 
ter of 2005, when the temperature did not reach the usual 
12.9°C minimum, which is in principle a recurrent charac- 
teristic of the BOUSSOLE site, because of the occurrence 
of deep mixing during winter. The latter is the mecha- 
nism for the formation of the deep waters in the western 
Mediterranean basin. The cooling of the surface waters 
started much earlier than usual in the fall of 2005 — as early 
as mid-December — instead of mid-January or February as 
in most other years. 

8.3 Phytoplankton Pigments 

Pigment samples for HPLC analysis have been collected 
at about 12 depths between the surface and 200 m at the 
BOUSSOLE site since July 2001. Surface samples at 5 and 
10 m depths were collected in triplicate. During certain 
surveys only surface samples were taken. In the follow- 
ing section, results from these analyses from July 2001 to 
December 2005 are briefly described. 

8.3.1 TChl a Seasonal Evolution 

During the winter and the start of the spring bloom, 
moderate to high values for [TChl a] have been observed, 
particularly for the month of February, when maximum 
values at the surface reached 1.34 in 2002 and 2.44 mg m -3 
in 2003 (Fig. 14). In February 2002, the water column ap- 
peared to have undergone strong vertical mixing just be- 
fore the start of the survey, with [TChl a] values averaging 
0.845 mg m -3 (±5%) down to 130 m. 

March and April are characterized by the spring bloom, 
which takes place in the near-surface waters, with [TChl a] 
values reaching approximately 5mgm -3 . In March 2003, 
a deep chlorophyll maximum (DCM) already began to 
form at a depth of 50 m. During summertime, surface 
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Fig. 13. A time series (July 2001 to June 2005) at the BOUSSOLE site of the potential temperature over the 
0-200 nr depth range (lighter to darker grays are associated with lower to higher temperatures, respectively). 
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Fig. 14. The variations in the determination of the total chlorophyll a concentration between July 2001 and 
June 2005 at the BOUSSOLE site. Two sampling depths are depicted: depths less than 5 nr are represented 
by the open circles, and depths greater than 5 nr and less than 10 nr are given by the solid circles. Each point 
is the average of all the measurements taken during a single cruise to the site. The vertical bars represent one 
standard deviation with respect to the average for each cruise. 

[TChl a] concentrations remain low (0.1-0.2nrgnr -3 and DCM at 50 nr tends to decrease towards the end of the 
as low as 0.05 mgm -3 in the summer of 2001). Figure 15, summer period. 

however, demonstrates that this trend is not reflected in Autumn marks the transition between stratified sunr- 
the deeper samples (below 10 nr) where important develop- nrer conditions and winter turbulence (mixing) of the water 
nrent of phytoplankton takes place during the oligotrophic column, with the progressive increase of surface [TChl a] 
conditions of this Mediterranean site. The thickness of the values. 
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Fig. 15. TChl a (solid circles) and divinyl Chi a (open circles) over the 0-200 nr depth range for the following 
days in 2003: a) 8 February, b) 17 March, c) 5 April, d) 27 June, e) 23 September, and f) 5 December. 
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Fig. 16. A time series of average abundances of pico- (circles), nano- (diamonds), and micro- (squares) 
phytoplankton for surface samples (less than 20 m) between July 2001 and October 2005. Each point is the 
average over all measurements taken during a cruise. The gray shading is arbitrarily assigned to enhance 
the interpretation of the time series. The light and middle gray portions show the dominant separation of 
the community composition (in terms of contribution to the TChl a concentration) into a smaller percentage 
of nricroplankton and a larger percentage of nanoplankton, respectively. The exceptions are denoted by 
the occasional reversals highlighted in the darkest gray. The seasonal contributions of the picoplankton are 
indicated by the bold line. 


8.3.2 Changes in Community Structure 

Six carotenoids and total chlorophyll b were used to par- 
tition the phytoplankton into three size classes based on 
the unique aspects of the pigments used to establish the 
classification scheme. The fractions of the Chi a concen- 
tration associated with each of the pico-, nano-, and micro- 
plankton classes (f p , f n , and f m , respectively) are derived 
using the following relationships (Uitz et al. 2006): 


f P 

fn 

and 


0.86 [Zea] + 1.01 [TChl b] 

‘ [5p] 

1.27 [Hex] + 0.35 [But] + 0.6[AUo] 


[DP] 


( 10 ) 

( 11 ) 


1.41 [Fuco] + 1.41 [Peri] 

; pp] ‘ ‘ 


( 12 ) 


dinoflagellates) and nanophytoplankton (characterized by 
prymnesiophytes and crysophytes) tend to vary in opposi- 
tion with picophytoplankton (characterized by cyanobac- 
teria at this site). The first are not a dominant class over 
the studied period, with minimum values in autumn. Nev- 
ertheless, at depth, microphytoplankton can reach impor- 
tant proportions relative to the other two classes, partic- 
ularly in spring. Nanophytoplankton is a dominant class 
in surface waters over most of the year, except during the 
autumn when picophytoplankton tend to take over. The 
DCM, however, is dominated by nanophytoplankton over 
the whole summer-autumn time period. 

The presence of divinyl Chi a, which is an indicator of 
prochlorophytes, has generally been observed during au- 
tumn at the DCM and even in surface waters during winter 
convection. In regards to the presence and concentrations 
of Chi a degradation products (e.g., chlorophyllide a and 
phaeophorbide a), they are most often below the detection 
limits of the HPLC method. 


where the diagnostic pigment, DP, is simply the sum of 
the seven pigment concentrations used to distinguish the 
three size classes. 

Figure 16 presents the variations of these three classes 
for near-surface (between 0-10 nr) waters between the mid- 
dle of 2001 and the end of 2005. Generally, for surface wa- 
ters, microphytoplankton (characterized by diatoms and 


8.4 AOPs from Ship Operations 

Figure 17 panels a, b, and c show the July 2001 to Au- 
gust 2004 time series of the reflectance ratio i?(A)/i?(560), 
for A=443, 490, and 510 nnr, respectively. Figure 18 shows, 
for the same time period, the diffuse attenuation coeffi- 
cients of the upper layers for the wavelengths 412, 490, 
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Fig. 17. A times series, from July 2001 to August 2004, showing three reflectance ratios: a) i?(443)/i?(560), 
b) i?(490)/i?(560), and c) i?(510)/i?(560). The open circles are the average of the in situ data (several points 
per cruises) over each cruise. The open diamonds are for the (HPLC) [TChl a] values, the solid circles are for 
the reflectance ratio as computed through a model that is based on the in situ [TChl a] values (Morel and 
Maritorena 2001), and the solid diamonds are the [TChl a] values that would be derived through an ocean 
color algorithm using the in situ reflectances ratio. 
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Fig. 18. A times series, from July 2001 to August 2004, showing three diffuse attenuation coefficients: a) 
Kd{ 412), b) A'd(490), and c) 1^(560). The open circles are the average of the in situ data (several points 
per cruises) over each cruise. The open diamonds are for the (HPLC) [TChl a ] values, and the solid circles 
are for the Kd values as computed through a model that is based on the in situ [TChl a] values (Morel and 
Maritorena 2001). 
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Fig. 19. A times series, from July 2001 to August 2004, showing three irradiance reflectances: a) i?(412), 
b) i?(490), and c) J?(560). The open circles are the average of the in situ data (several points per cruises) 
over each cruise. The open diamonds are for the (HPLC) [TChl a] values, and the solid circles are for the 
R values as computed through a model that is based on the in situ [TChl a] values (Morel and Maritorena 
2001 ). 
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and 560 nrn. Finally, Figure 19 shows the reflectances, 
R = E u (0~) / Ed(0~), for the same wavelengths. For all 
of these figures, the open circles are for the in situ data, 
the open diamonds are for the total chlorophyll a concen- 
tration determined from HPLC analysis, the solid circles 
are for the AOPs (either R , the ratio of R at two wave- 
lengths, or Kd) as computed through a model that uses 
the in situ Chi a concentration as the input variable, and 
the solid diamonds are the Chi a concentration that would 
be derived through satellite algorithms using the in situ 
reflectance ratios. 

The difference between the open and solid circles, there- 
fore, represents the anomaly of the measured AOPs as 
compared to what is predicted by standard (global) bio- 
optical models, considering the in situ Chi a concentra- 
tion. Similarly, the difference between the open and solid 
diamonds represents the uncertainty in the Chi a concen- 
tration that is obtained when satellite algorithms, which 
were derived from global data sets and are based on average 
global relationships between the reflectance ratios and the 
Chi a concentration, are applied to waters with anomalous 
optical properties. 

One striking feature is the large anomaly of the blue- 
to-green ratio in the summer of 2001, with lower than ex- 
pected blue-green ratios, as already identified and quan- 
tified in Claustre et al. (2002) for the 1999 Productivity 
des Systemes Oceaniques Pelagiques (PROSOPEf) cruise. 
The reason for the anomaly was attributed to the pres- 
ence of desert dust particles in the water, which have an 
enhanced absorption in the blue and an enhanced scat- 
tering in the green parts of the spectrum (Claustre et al. 
2002). This anomaly is, however, less important during 
the summer of 2002, it nearly vanishes during the summer 
of 2003, and appears again in the summer of 2004. The 
transient character of the anomaly has yet to be confirmed 
and explained. 

Examination of the relationships for the diffuse attenu- 
ation coefficients, Kd{ A), in Fig. 18 reveals that their values 
in the summer of 2001 are largely above the modeled val- 
ues across the entire spectrum, at least from 412-560 nm, 
and that this difference is disappearing in summer of 2003. 
Because Kd is largely determined by the absorption prop- 
erties of the medium, these observations mean that an ex- 
cess absorption exists at all wavelengths in the July to 
September period in 2001. 

The same feature does not appear in the reflectances 
(Fig. 19), which are lower than the modeled values in the 
blue, again during the summer of 2001, yet greater than 
the modeled values in the green. Because R is, to the 
first order, proportional to backscattering and inversely 
proportional to absorption, these observations confirm the 
presence of an excess of backscattering in the green part of 


f Translated as productivity of pelagic oceanic systems. De- 
tailed information about the PROSOPE cruise is available 
from the http//www . obs-vlf r . fr Web site. 


the spectrum. Again, the transient nature of this anomaly 
remains to be confirmed and explained. 

Another anomaly is seen in February 2003 with larger- 
than-expected blue-green ratios, which is in opposition to 
the summer anomaly. During the February 2003 cruise, the 
Chi a concentration was as large as 2 mg nr -3 . The reason 
for this behavior might be in the low scattering coefficients 
that are typical of a fresh phytoplankton bloom, where the 
proportion of large healthy cells characterized by a low 
backscattering efficiency is high, and the contribution of 
small detritus particles is low. 

Figure 20 shows the seasonal and interannual evolution 
of Kd for wavelengths 412, 490, and 560 nm. The time se- 
ries was established from three years of measurements car- 
ried out during the (approximately) monthly maintenance 
cruises. This parameter is largely determined by the ab- 
sorption coefficient, and it is thus a good indicator of the 
concentration of phytoplankton and associated colored dis- 
solved substances. The seasonal cycle appears very clearly, 
in particular in 2003 and 2004: 

• It starts with an intense mixing of water in winter 
(maximum in February) resulting in almost homo- 
geneous water properties and weak Kd values; 

• It continues with the start of the spring bloom in 
the March- April time period; and 

• It then evolves gradually to oligotrophic conditions 
(with a minimum in Chi a concentrations around 
September), which is characterized by a peak in at- 
tenuation at about 50 nr depth (coincident with the 
DCM). 

For Chi a concentrations lower than 0.1 mg nr -3 at the sur- 
face during the latter oligotrophic period, the global sta- 
tistical relationships (Morel and Berthon 1989) predict a 
DCM at approximately 100 nr. Its development at a lesser 
depth at the BOUSSOLE site is due to the fact that the 
site is in the center of a cyclonic circulation, which is 
dynamically represented as a donre-like structure that is 
maintained by significant vertical advection (the vertical 
ascent of deeper water transports the DCM to a shallower 
depth) . 

The interannual variability observed here for Kd seenrs 
directly related to the changes of the physical framework. 
Indeed, the depth of the DCM seenrs to increase by 2001 
to 2004, at the same time as the winter mixing seenrs to 
intensify and the spring bloom being more intense. The 
oligotrophic layer of water characterized by low values of 
Kd in summer is deeper in 2003 than in 2001 and 2002. 
An intense mixing in winter on the one hand brings more 
nutrients to the surface, supporting a stronger bloom, and, 
on the other hand, more efficiently eliminates the colored 
dissolved substances accumulated at the surface during 
summer. This phenomenon has already been shown for 
dissolved organic carbon (Copin- Montegut and Avril 1993 
and Avril 2002) and could explain the changes observed 
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Fig. 20. A time series of the vertical distribution of the diffuse attenuation coefficient ( Kd ) derived from 
SPMR Ed vertical profiles at a) 412 nm, b) 490 nrn, and c) 560 nm. The small black bars protruding upwards 
from the rr-axis and into the contour maps denote the time of the ship sampling. The units for Kd are per 
meter. 
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Fig. 21. A time series of data collected in 2003 of the absorption coefficients at 440 nrn for total particles 
(open squares), phytoplankton (open circles) and detritus (open diamonds). The absorption coefficient of 
phytoplankton at 676 nrn is also depicted (solid circles). Two different scales for the y-axis are shown: a) a 
linear scale, and b) a logarithmic scale in order to provide a better discrimination of the low values during 
the summer months. 


in the so-called anomalies of the marine reflectance as dis- 
cussed above. 

The interannual variability of the spring bloom, how- 
ever, is somewhat distorted by an insufficient sampling in 
spring 2002 and 2003 (bad weather). Consequently, it will 
probably be necessary to rebuild part of the record by using 
the vertical profiles of the pigments and the values derived 
from the satellites, coupled with the relationships between 
Chi a and K d that will be established from the existing 
data. 

8.5 IOPs from Ship Operations 

An annual cycle of the particulate absorption coeffi- 
cients at 440 nrn is presented in Fig. 21, where the differ- 
ent curves are for total particles (a p , directly measured on 
filtered samples), and both the phytoplankton and detrital 
particles ( a $ and ad, respectively). The latter are derived 
through a numerical decomposition scheme (Bricaud and 
Stramski 1990) as presented in Sect. 7.2. An important 
feature here is the high contribution of the absorption by 
detrital particles, at least 50% of the total particulate ab- 
sorption, from about the end of June to the beginning of 
September. 

Figure 22 shows examples of the absorption spectra 
by particles, either for the total particles or for the phy- 
toplankton or detritus only. The dominant role of phyto- 
plankton in the absorption budget during the spring bloom 
clearly appears, as well as the progressive increase of that 
part of absorption that is due to detritus in summer. The 


exact cause of the background absorption in the ultravi- 
olet (UV) that clearly appears when absorption decreases 
at other wavelengths (see the summer spectrum) remains 
to be identified. 

Figure 23 presents examples of the total attenuation, 
scattering, and absorption coefficients for the surface (0- 
10 m) layer (from the AC-9+). Note the clear identification 
of the absorption features on the scattering spectrum, and 
the change in its slope from spring to summer (the increas- 
ing slope in summer corresponds to the large contribution 
of detritus as compared to phytoplankton alone). 

Selected observations of the total scattering coefficient 
at 550 nm, as determined from AC _ 9+ measurements, are 
displayed in Fig. 24 as a function of the TChl a concen- 
tration, along with the values derived from the Loisel and 
Morel (1998) relationship. What can be noticed is that the 
organization of the in situ data with respect to the model is 
providing some clues for their interpretation. For instance, 
the three lowest values, which correspond to summer, are 
above the model line; this may be linked to the increasing 
proportion of detritus during this season. In contrast, the 
largest value, which corresponds to the spring bloom, is 
largely below the model line; this may indicate the domi- 
nance of large healthy cells with low scattering efficiencies. 
The other points in the graph, which are closer to the 
model line, correspond to autumn. This example and the 
tentative interpretation of the data distribution are pro- 
vided to illustrate the type of investigations that will be 
performed with the full data set that is being collected for 
the BOUSSOLE activity. 
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Fig. 22. Spectra of the absorption coefficient of particies at 5 m for 2003: a) 8 February, b) 17 March, c) 29 
May, and d) 13 July. The [TChl a] values are 1.85, 0.52, 0.15, and 0.05 mg m~ 3 , respectively. 
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Fig. 23. A09+ spectra of the total attenuation (circles), absorption (diamonds), and scattering (squares) 
coefficients for the near-surface layer (10 m) for four different time periods: a) 7 February 2003, b) 5 April 
2003, c) 27 May 2003, and d) 31 July 2004. 
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Fig. 24. The total scattering coefficient at 550 nrn 
as a function of the TChl a concentration. Solid 
circles correspond to in situ values ( AC-9+) , and the 
solid and dashed lines are 6(550) values computed 
from Loisel and Morel (1998) and a factor of two 
around the modeled values, respectively. 

8.6 AOPs from Buoy Data 

Before discussing the time series of the reflectances de- 
rived from the buoy data, a comparison is performed with 
the reflectances derived from the SPMR. Although a lim- 
ited data set has been used here (not the full time series 
of buoy and SPMR data), the comparison (Fig. 25) shows 
a good agreement, with a small bias of approximately 2%. 



[Lw] n Buoy [|iW cm -2 nnrr 1 sr 1 ] 

Fig. 25. A comparison of normalized water-leaving 
radiances from the SPMR and buoy for the wave- 
lengths indicated. 


A least squares fit to all the Fig. 25 data (dashed line) 
has a slope of 0.93, a y - intercept of 0.01, and a coefficient 
of determination of r 2 = 0.98. Further work is needed to 
unambiguously attribute this bias to either the SPMR or 
the buoy. 

Figure 26 presents a one-month time series of the re- 
flectance at 443 nm and 560 nrn, corresponding to the de- 
velopment of the spring bloom in 2004 (March 4 to April 7) . 
The TChl a concentration ranges from about 0.2 mg nr -3 
at the beginning of the sequence to about 3.5 mg nr -3 at 
the end, and the reflectance in the blue and green domains 
change accordingly. The reflectance goes from about 3% to 
less than 1% at 443 nm, and from about 1-2% at 560 nm. 
Large daily changes also occur, which result from the com- 
bined effects of changes in the water optical properties and 
the daily solar cycle. The latter can be easily removed 
through an appropriate normalization (Morel and Gentili 
1991). 

Figure 27 shows a longer time series of the reflectance 
in the blue (443 nm) , which actually corresponds to all the 
data collected by the buoy in 2005. The full range of vari- 
ability spans about 1-5% in reflectance. The depression 
caused by the spring phytoplankton bloom is maximum 
around mid- April, and the maximum of reflectance occurs 
at the end of July, when the TChl a concentration is close 
to its annual minimum (about 0.1 mg m -3 ). Note that the 
reflectance significantly decreases during August, whereas 
the TChl a concentration does not increase and remains 
close to O.lmgm -3 . This phenomenon is another illus- 
tration of what has been already discussed in Sect. 8.4 
( anomalies in the blue-to-green ratio). 

8.7 IOPs from Buoy Data 

Figure 28 displays the record of the phytoplankton fluo- 
rescence (dashed curve) and of the beam attenuation coef- 
ficient at 660 nm (solid curve), during a one-month period 
corresponding to the development of the spring bloom in 
2004 (March 4 to April 7; same as in Fig. 26). The TChl a 
concentration is from about 0.2 mg m -3 at the beginning 
of the sequence to about 3.5 mg m -3 at the end. 

The fluorescence varies in a complex manner. For in- 
stance, the depression of the fluorescence signal around 
noon, which depends on the actual irradiance at the sea 
surface, is well pronounced. It may actually lead to vari- 
ations within a day that are larger than the day-to-day 
change of the average fluorescence. A practical conse- 
quence is that the signal during the day cannot be used 
as a proxy of the Chi a concentration. Night values are 
better adapted for that purpose, and indeed they seem to 
be well correlated to the evolution of the Chi a concentra- 
tion (not shown here). 

The beam attenuation coefficient at 660 nm is a proxy 
of the particle load, and varies closely with the fluores- 
cence values recorded at night. The bloom produces a 
ten-fold increase in c(660), from 0. 2-2.0 nr -1 , mostly oc- 
curring between days 90 and 95. Daily changes are also 
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Fig. 26. An approximately 4- week buoy time series of R( 443) and i?(560), indicated as black and gray circles, 
respectively. The thin vertical lines indicate solar noon. 
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Fig. 27. A time series of all the reflectance data in the blue part of the spectrum (443 nrn) collected by 
the buoy around noontime in 2005. The black squares at the top of the plot denote the occurrences of the 
monthly maintenance cruises. 
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Fig. 28. An approximately 4- week buoy time series of the beam attenuation coefficient at 660 nrn (black 
circles) and the phytoplankton fluorescence (gray circles). The latter is presented in relative units. 


clearly shown, with a great variety of shapes, however, and 
sometimes with very fast increases (day 95 for instance) or 
dramatic declines (day 96). The maximum is usually oc- 
curring around noon or a few hours after noon. Inherently 
to Eulerian measurements at a fixed moored station, these 
changes are the result of the local phytoplankton growth 
and of lateral advection of water masses. 

8.8 Sun Photometer Data 

Figure 29 illustrates the overall data set that has been 
collected up to now at the Cap Ferrat sun photometer 
site. What is shown is the aerosol optical thickness as de- 
rived by the AERONET project through their standard 
procedures. The monthly AOT averages are usually min- 
imum in January or February, and several relative max- 
ima occur in spring (March or April), summer (usually 
July), and fall (October-November). The corresponding 
monthly averages of the Angstrom exponent (not shown) 
are always above 1.0m -1 , which indicates the persistence 
of small particles throughout the year. 

The ability of the inversion method, which was devel- 
oped as part of the BOUSSOLE project and uses ground- 
based data to get aerosol optical properties (Sect. 7.7), is 
demonstrated using days corresponding to MERIS over- 
passes. The results are presented in Fig. 30 for sequences 
of principal planes corresponding to high solar zenith an- 
gles (typically above 70°) to better cover the backscatter- 
ing region of measurements and to include the scattering 
angle of 150°. Table 8 presents the values of the retrieved 


refractive index for these sequences. The agreement be- 
tween measurements and predictions is good, especially in 
the backward direction for sky radiances. The degree of 
polarization is also accurately retrieved, especially at 90°, 
which is of interest for classifying particle types. 

These examples demonstrate it is possible to obtain 
some information on the refractive indices of aerosols by 
using the inversion method and, thus, to get the aerosol 
phase function. Such information can be used to predict 
the angular distribution of the radiance and the degree of 
polarization, thus, it is likely that ground-based photome- 
ter measurements can be used for the vicarious calibration 
of satellite ocean color sensors (Sect. 10). 

9. SATELLITE VALIDATION 

Match-up analyses of the water-leaving radiances, as 
well as for the Chi a concentrations, have been performed 
for the MERIS, SeaWiFS, and MODIS-A data products. 
The analyses all follow the same procedure, wherein a 5x5 
pixel box is extracted from the level -2 product (see below), 
and directly compared to the same quantity as derived 
from the in situ data (ship or buoy). For the SeaWiFS and 
MODIS~A sensors, the criteria used for a given comparison 
point to be excluded in the final analysis are described by 
Bailey et al. (2000). They are essentially the same for 
MERIS data, except the flags are not identical. 

The combination of data availability, orbit character- 
istics, and the quality criteria used in match-up process 
yields 68 matchups for MERIS (data from September 2003 
to December 2005), 89 matchups for SeaWiFS (data from 
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Fig. 29. The full time series of the AOT at the AERONET Cap Ferrat coastal site with monthly averages 
shown as the larger open circles. The corresponding calendar years are given along the top of the plot. 
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Fig. 30. Three examples of reconstructing sky radiances (in units of microwatts per square centimeter per 
nanometer per steradian) in the principal plane and at three wavelengths (curves), as compared to the direct 
measurements of the sun photometer (symbols), from which the inversion procedure (Sect. 7.7) infer an aerosol 
type: a) 11 July, a) 7 September, and c) 26 September 2002. The corresponding polarization rate at 870 nm 
is shown in panels d— f. 
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Table 8. The aerosol optical properties returned by the inversion algorithm for three days of measurements. 
These parameters were used in the computations presented in Fig. 30. The complex refractive index, m r — jrrn, 
aerosol optical depth at 675 nrn, and the Angstrom exponent, a, measured by the AERONET radiometer are 
presented along with the solar zenith angle, 6 S . 


Year 

Date 

Time [GMT] 

m r 

- jmi 

r a (675) 

a 

0s [°] 

2002 

11 July 

1736 

1.44 - 

jO.0049 

0.0955 

1.562 

72.34 

2002 

7 September 

1627 

1.41 - 

jO.OOOOl 

0.1143 

1.686 

72.79 

2002 

26 September 

1567 

1.51 - 

j 0.0090 

0.0385 

1.459 

72.82 


September 2003 to December 2004) and 65 matchups for 
MODIS-A (data from September 2003 to May 2005). The 
data used are a function of the satellite mission. For 
MERIS, level-2 1 km reduced resolution data processed by 
the MERIS Ground Segment Simulator (MEGS) prototype 
version 7.4.1 are used. For SeaWiFS, level -2 merged lo- 
cal area coverage (MLAC) data from reprocessing 5 (com- 
pleted March 2005) are used, and for MODIS-A, level-2 
global area coverage (GAC) data from reprocessing 1 (com- 
pleted February 2005) are used. 

The MERIS data provide the water-leaving radiance 
reflectance, p w , defined as (again ignoring the spectral de- 
pendence for brevity): 


Pw(e v ,9 a ,A<i>) 


Lw(9 v ,0 s , A<p) 

E s (e s ) 


(13) 


where neither Lw nor Ed are normalized to a zenith Sun 
or to a nadir viewing angle. It is worth noting as well that 
the MERIS data have still not been vicariously calibrated. 

The SeaWiFS and MODIS-A data provide the fully 
normalized water-leaving radiance, [Lw(A)] which are 
defined (without the spectral dependence) as: 


[L 


Lw 


% 


W Ur - 


t d (d s )cose s 


f n(//<2), (14) 


where JJ is a geometrical factor accounting for all of the 
refraction and reflection effects at the air-sea interface, 
Jfo = $1(0' = 0), and fn (f/Q) provides the functional de- 
pendence on the bidirectional nature of the light held. The 
latter is defined as 

-l 

, ( 15 ) 

where C a is the Chi a concentration, W is wind speed, Q is 
the ratio of upward irradiance to upwelling radiance (sub- 
script 0 means for a Sun at zenith and a nadir view) , / is a 
factor relating the irradiance reflectance R to the inherent 
optical properties, and other symbols are describing the 
geometry of the problem. The transformations needed to 
derive fully compatible quantities between the satellite and 
the in situ measurements (i.e. , accounting for the bidirec- 
tional nature of the light held) are not detailed here. 


M//Q) = 


fo(W,C a ) 
Qo(W,C a 


f(e 8 ,W,C a ) 
Q(0 a ,0',A<t>,W,C a ) 


Concerning the water-leaving radiances or rehectances, 
results are displayed as scatterplots of the satellite-derived 
versus in situ measurements in Fig. 31 for MERIS, Sea- 
WiFS, and MODIS-A (all wavelengths presented together). 
The associated statistics are provided separately for each 
wavelength, as well as for the f ull data set, in Tables 9-11, 
for MERIS, SeaWiFS, and MODIS-A, respectively. 

The overall agreement between the satellite-derived and 
the in situ rehectances or radiances is similar for the three 
satellites, and for SeaWiFS and MODIS-A, it is close to 
what has been already reported elsewhere. Consequently, 
a table-by-table discussion of the statistics within the ta- 
bles is not presented, but some specihc points are noted: 

• The MERIS data products for the red bands (620, 
665, and 681 nrn) are overestimated by a factor of 
about two, which is not seen for the other satellite 
data products; 

• The MERIS data products for the blue wavelengths 
appear to be overestimated; 

• The MODIS-A data products exhibit the smallest 
amount of dispersion around the 1:1 line; and 

• The SeaWiFS and MODIS-A data products exhibit 
a slight underestimation in the blue domain (the 
slopes of the overall fit are slightly less than 1). 

The TChl a matchups for the three sensors are pre- 
sented in Fig. 32, where all available comparison points 
are shown regardless of the environmental conditions, and 
in Fig. 33, where only a quality-assured subset of data 
is shown based on the criteria used for the reflectance 
matchups. The number of match-up points in Figures 32 
and 33, respectively, are 34 and 9 for MERIS, 52 and 14 
for SeaWiFS, and 39 and 9 for MODIS-A. 

10. VICARIOUS CALIBRATION 

The TOA vicarious radiometric calibration is a chal- 
lenging process, for which high-quality data from several 
sources have to be brought together, to form input to, 
and boundary conditions of, radiative transfer calculations. 
These calculations are supposed to provide the TOA ra- 
diances, which are compared to the same quantities as 
measured by the spaceborne sensor. The principle of this 
calibration is given below, before a tentative performance 
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0 0.01 0.02 0.03 0.0001 0.001 0.01 0.03 

In Situ p w { A) In Situ p w ( A) 




In Situ [L w (A)] N In Situ [L w (A)] N 



0 0.4 0.8 1.2 1.6 0.003 0.01 0.1 1 3 


In Situ [L w (A)] N In Situ [L w (A)] N 

Fig. 31. Satellite matchups at the BOUSSOLE site for a) MERIS p w ( A) values (412, 443, 490, 510, 560, 
670, and 683 nm); b) SeaWiFS [Lw(A)] n , values (412, 443, 490, 510, 555, and 670nm); and c) MODIS~A 
[Zw(A)] n values (412, 443, 488, 510, 551, 670, and 683 nm). The solid line is the 1:1 line. Logarithmic scales 
for panels a c are shown in panels d— f, respectively, in order to magnify the low values in the red domain. 
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Table 9. The match-up statistics for MERIS. The statistical parameters include the relative percent difference 
(RPD) plus four linear fit parameters: the coefficient of determination, r , the slope, m; the y-intercept, y 0 ; and 
the root mean square (rms) error. 


A 

[nrn] 

No. 

Pts. 

Mean 

Ratio 

RPD 

[%] 

Linear Fit Parameters 
r 2 m y 0 rms 

In Situ 

Min. Max. 

Satellite 
Min. Max. 

412 

29 

1.26 

25.69 

0.638 

0.84 

0.0055 

0.00468 

0.00512 

0.02420 

0.00760 

0.02784 

443 

55 

1.10 

10.42 

0.722 

0.86 

0.0033 

0.00275 

0.00415 

0.02357 

0.00615 

0.02372 

490 

55 

1.09 

8.62 

0.750 

0.79 

0.0037 

0.00181 

0.00448 

0.01772 

0.00621 

0.01852 

510 

55 

1.13 

12.58 

0.437 

0.56 

0.0053 

0.00169 

0.00496 

0.01277 

0.00653 

0.01281 

560 

55 

1.14 

14.36 

0.237 

0.44 

0.0036 

0.00111 

0.00332 

0.00815 

0.00401 

0.00869 

620 

7 

1.93 

92.53 

0.074 

0.76 

0.0008 

0.00076 

0.00051 

0.00094 

0.00070 

0.00213 

665 

55 

3.42 

241.51 

0.022 

0.46 

0.0008 

0.00074 

0.00013 

0.00076 

0.00029 

0.00186 

681 

55 

5.32 

431.90 

0.509 

1.24 

0.0007 

0.00083 

0.00001 

0.00163 

0.00033 

0.00261 

709 

7 

2.00 

99.62 

0.228 

1.24 

0.0002 

0.00042 

0.00013 

0.00061 

0.00007 

0.00111 

All 

373 


0.938 

1.02 

0.0009 




Table 10. The match-up statistics for SeaWiFS. 


A 

[nrn] 

No. 

Pts. 

Mean 

Ratio 

RPD 

[%] 

Linear Fit Parameters 
r 2 m y 0 rms 

In Situ 
Min. Max. 

Satellite 
Min. Max. 

412 

48 

0.95 

-5.10 

0.770 

1.16 

-0.20 

0.169 

0.367 

1.545 

0.263 

1.548 

443 

98 

0.99 

-0.96 

0.795 

0.96 

0.03 

0.119 

0.280 

1.430 

0.314 

1.525 

490 

98 

1.00 

-0.20 

0.673 

0.78 

0.19 

0.100 

0.315 

1.194 

0.384 

1.189 

510 

98 

0.99 

-1.49 

0.403 

0.57 

0.25 

0.078 

0.335 

0.848 

0.359 

0.796 

555 

98 

1.01 

1.08 

0.123 

0.37 

0.19 

0.054 

0.187 

0.498 

0.173 

0.436 

670 

95 

1.37 

36.70 

0.026 

0.22 

0.02 

0.014 

0.007 

0.065 

0.004 

0.061 

All 

535 


0.950 

0.97 

0.01 




Table 11. The match-up statistics for MODIS-A. 


A 

[nrn] 

No. 

Pts. 

Mean 

Ratio 

RPD 

[%] 

Linear Fit Parameters 
r 2 m y 0 rms 

In Situ 
Min. Max. 

Satellite 
Min. Max. 

412 

42 

0.95 

-4.62 

0.864 

1.03 

-0.06 

0.120 

0.405 

1.605 

0.282 

1.590 

443 

85 

0.98 

-2.13 

0.824 

0.91 

0.07 

0.094 

0.310 

1.495 

0.286 

1.585 

488 

85 

1.00 

-0.10 

0.795 

0.79 

0.19 

0.069 

0.354 

1.174 

0.372 

1.228 

551 

85 

1.08 

8.31 

0.242 

0.59 

0.15 

0.049 

0.196 

0.380 

0.208 

0.432 

667 

84 

1.48 

48.44 

0.003 

0.06 

0.02 

0.010 

0.005 

0.048 

0.005 

0.040 

All 

381 


0.975 

0.96 

0.02 




budget is proposed and preliminary results are shown for 
the MERIS sensor. 

10.1 Principles 

Two main vicarious f calibration paths exist to produce 
ocean color products of the desired accuracy, i.e. , water- 
leaving radiances within an uncertainty of about 5% in the 
blue for an oligotrophic ocean (Gordon 1997, and Antoine 


f The true , direct, calibration in principle consists in measur- 
ing the signal from a well-known standard, and is performed 
before launch for spaceborne sensors. After launch, only in- 
direct ( vicarious ) procedures are possible. 


and Morel 1999). The first one is usually referred to as 
vicarious calibration, and consists in forcing the satellite- 
derived water-leaving radiances to agree with a set of in 
situ water- leaving radiances (match-up analyses) . A set of 
vicarious calibration coefficients is obtained, which is ap- 
plied to the TOA radiances measured by the sensor. The 
second procedure, which is also an indirect (vicarious) cal- 
ibration is sometimes referred to as a vicarious radiometric 
calibration, and consists in simulating the TOA signal that 
the sensor should measure under certain conditions, and to 
compare it to the measured signal. 

One of the difficulties of the first type of vicarious cal- 
ibration is that it is dependent upon the procedure used 
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In Situ TChl a [mg rrr 3 ] In Situ TChl a [mg nrr 3 ] In Situ TChl a [mg rrr 3 ] 


Fig. 32. In situ (HPLC) TChl a concentration matchups with respect to satellite-derived values at the 
BOUSSOLE site for a) MERIS (N = 34), b) SeaWiFS (N = 52), and c) MODIS-A (N = 39). The temporal 
coverage is from July 2001 to the end of 2005 for MERIS and MODIS-A, and 2004 for SeaWiFS. 



In Situ TChl a [mg rrr 3 ] In Situ TChl a [mg nr 3 ] In Situ TChl a [mg nr 3 ] 


Fig. 33. The quality- assured in situ (HPLC) TChl a concentration matchups shown in Fig. 32 for a) MERIS 
(N = 9), b) SeaWiFS (N = 14), and c) MODIS-A (N » 9). 


for the atmospheric correction of the TOA observations. 
Even if it is admittedly less dependent upon the selected 
set of in situ water-leaving radiances, these measurements 
also contribute to the final accuracy. Besides the fact that 
atmospheric measurements are not needed, the advantage 
of this technique is that the marine signals delivered by 
several sensors using different atmospheric correction al- 
gorithms can be cross-calibrated! provided the same set 
of in situ water-leaving radiances is used to perform the 
vicarious calibration (presently true for the SeaWiFS and 
MODIS-A sensors). 

An inconvenience of the vicarious radiometric calibra- 
tion is that it requires a set of in situ measurements which 
are usually difficult to collect, among other things, because 


f Except in situations where the atmosphere would be signifi- 
cantly different from the atmosphere present when perform- 
ing the in situ measurements used for the vicarious calibra- 
tion exercise. 


a high accuracy is required. In addition to the in- water 
measurements of the water-leaving radiances, this data set 
includes sea state and atmospheric pressure, ozone concen- 
tration, aerosol optical thickness, aerosol type, and even 
aerosol vertical profile if the aerosols reveal to be absorb- 
ing. If this data set is successfully assembled, the ad- 
vantage of the vicarious radiometric calibration is that it 
is independent of the atmospheric correction algorithms, 
so that the TOA signals of various sensors can be cross- 
calibrated. Then it is up to any user to apply its preferred 
atmospheric correction to these TOA signals. The marine 
signals in that case might be inconsistent if significant dif- 
ferences exist in the various atmospheric corrections. 

The greatest difficulty of the vicarious radiometric cal- 
ibration lies in the estimation of the aerosol optical thick- 
ness, phase function, and single scattering albedo. These 
parameters are accessible through the inversion of sun pho- 
tometer measurements, yet uncertainties inevitably occur 
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when applying such methods, for example because of mul- 
tiple scattering, perturbations from the ground reflectance, 
or uncertainties in the photometer calibration. Assem- 
bling all data needed for these vicarious calibration experi- 
ments are often compromised only because aerosol param- 
eters are not accurate enough. The principle is illustrated 
on Fig. 34. A practical example is provided later on in 
Sects. 10.3 and 10.4. 


• Atmospheric Pressure 

• Sea State (Wind Speed) 

• Aerosol Optical Thickness 

• Aerosol Type P(8)\ 

• Ozone Concentration 

• Geometry 



Lt 

Measured 


L 


path 


TLg T 


© 


+ td,Lw 



©Selection of the relevant geometry and 
sea state so that: 

TLg — tdL wc — 0 

©Measured with in-water radiometers or 
assumed to be zero in the near infrared 


Fig. 34. The overall logic of the TOA vicarious ra- 
diometric calibration process. L t is the total (TOA) 
radiance, L g is the radiance contribution from sun 
glint, and L wc is the radiance from white caps. The 
direct and diffuse atmospheric transmittances are 
denoted by T and td , respectively. 


a given process to the error budget was missing or insuffi- 
ciently known, this contribution has been assigned a some- 
what optimistic value. 

Each term that may impact the modeling of the TOA 
radiance ( L t ), which is the sum of the atmospheric path 
radiance (L pa th) and of the product of the water-leaving 
radiance (Lw) by the atmospheric diffuse transmittance 
(td), is examined below. At the end, the uncertainties on 
all terms are considered as random, gaussian distributed, 
independent one of each other, so that the final error bud- 
get is computed as the square root of the sum of the squares 
of the individual error terms. The extent to which devi- 
ation from this hypothesis may impair the correctness of 
the final error budget has not been assessed. 

The various uncertainty estimates discussed below, as 
well as the overall uncertainty budgets, are summarized in 
Tables 12 and 13. 

Table 12. Theoretical and tentative uncertainty 
budget for vicarious calibration in the near infrared. 

In this case, Lw is assumed to be 0 (i.e. , the wave- 
length is 865 nrn in the near infrared). 


Source of Uncertainty 

U L t 

AOT 

2.00 

Aerosol Type 

5.00 

Atmospheric Pressure 

0.20 

Sea Surface State 

0.50 

Ozone 

0.00 

Radiative Transfer 

0.00 

Overall Uncertainty 

5.41 


The Rayleigh scattering above clear waters method ac- 
tually follows the same principle as the full radiometric vi- 
carious calibration (the scheme displayed just above), but 
by using assumptions about, instead of measurements of, 
the various input parameters. This is particularly critical 
for the water-leaving radiances and the aerosol parameters. 
The final accuracy of this method cannot be better than 
that of the vicarious calibration procedure, except in a sit- 
uation where the measurements used in the latter would 
be of very poor quality (or even erroneous). 

10.2 Performance Budgets 

It is in practice impossible to establish an a priori full 
performance budget, which should account for both un- 
predicted or unknown experimental sources of error (such 
as unidentified calibration deficiencies, difficulties encoun- 
tered at sea, etc.) and theoretical limitations of the mea- 
surements themselves. The former are not totally known 
before starting the measurement sessions. 

The budget presented below is, therefore, a tentative 
one, and any number provided in terms of percent error 
should be considered with caution. When the informa- 
tion necessary to confidently estimate the contribution of 


10.2.1 In- Water Measurements 

The routine operations on the BOUSSOLE buoy should 
provide the measurements of Ed, E u , and L u (nadir) at 5 
and 9 m. The basic operation consists of estimating the 
water- leaving radiance, Lw, from these in- water measure- 
ments, and in particular from L u . Evaluation of above- 
water techniques for the measurement of Lw will also be 
tentatively set up (using a SIMBADA instrument), and 
will be examined later (Sect. 10.2.2). The uncertainty in 
L u is estimated to be about 5% when carefully performing 
the measurements with well-calibrated instruments. 

10.2.1.1 Extrapolation to Null Depth 

Here the problem is to compute the diffuse attenuation 
coefficient for the radiance along the nadir direction, A©, 
from the measurements performed at two depths on the 
buoy (nominally 5 and 9 m). Denoting this value as K ^ 9 , 
it is in principle used to compute the upwelling nadir ra- 
diance just below the surface, L u ( 0"), through 

L u ( 0-) = L u (Z)e K i 9{z \ (16) 

where Z is about 5 m. 


50 




Antoine et al. 


Table 13. Theoretical and tentative uncertainty budgets for in-water (left column), above-water (middle 
column), and Rayleigh scattering above clear water vicarious calibration (right column) experiments. For all 
three activities, Lw is 15% of L t , u L is the uncertainty expressed as a percentage of Ly y, and u L is the 
uncertainty expressed as a percentage of L t . 


In-Water 

U r 

L w 

U T 

L t 

Above-Water 

U T 

L w 

U T 

L t 

Rayleigh Scattering 

U r w 

U T 

L t 

Uncertainty Source 

[%] 

[%] 

Uncertainty Source 

[%] 

[%] 

Uncertainty Source 

[%] 

[%] 

Measurement of L u 

5.0 

0.75 

Measurement of Lw 

10.0 

1.50 

Assumption about Lw t 

30.0 

4.50 

Extrapolation to 0 + 

3.0 

0.45 







Bidirectionality 

2.0 

0.30 

Bidirectionality 

2.0 

0.30 




Atmos. Transmittance 

2.0 

0.30 

Atmos. Transmittance 

2.0 

0.30 

Atmos. Transmittance 

3.0 

0.45 

AOT 


1.00 

AOT 


1.00 

AOT Extrap, from NIR 


2.00 

Aerosol Type 


2.00 

Aerosol Type 


2.00 

Aerosol Type 


5.00 

Atmospheric Pressure 


0.20 

Atmospheric Pressure 


0.20 

Atmospheric Pressure 


0.50 

Sea Surface State 


0.50 

Sea Surface State 


0.50 

Sea Surface State 


1.00 

Ozone 


0.50 

Ozone 


0.50 

Ozone 


0.50 

Radiative Transfer 


0.00 

Radiative Transfer 


0.00 

Radiative Transfer 


0.00 

Overall Uncertainty 


2.55 

Overall Uncertainty 


2.82 

Overall Uncertainty 


7.14 


f With a priori knowledge of the area under examination. 


Three important points intervene here: a) the inter- 
calibration of the instruments at the two depths, which 
is assumed to be properly established; b) the correct es- 
timation of the depth of the measurements'!; and c) the 
adequacy of K'f 9 to perform the extrapolation of L u from 
5 m to just below the sea surface. 

For clear waters (i.e. , Kl about 0.02 m -1 in the blue), a 
large uncertainty of 1 m on the exact depth of measurement 
would lead to an error equal to e h ' L , i.e., a 2% uncertainty 
on the estimation of L u ( 0”). Conversely, if it is assumed 
that Z is correctly estimated, an uncertainty of 10% on Kp 
would lead to an uncertainty equal to e -°- 1K z.z , p e ^ an 
uncertainty of 1% on L u ( 0") if Z is taken equal to 5m (or 
2% if Z is 10 m). These numbers become 2.5% and 5% for 
K l = 0.05, i.e., for mesotrophic waters, also typical of the 
BOUSSOLE site when the TChl a concentration reaches 
about 0.3 mgm -3 . 

The attenuation coefficient for the upwelling nadir ra- 
diance is known to vary a bit within the very upper lay- 
ers of the ocean (2-5 m), so that it is timely to examine 
whether or not remains the relevant coefficient to 

perform the extrapolation of L u measured at 5 m to the 
0“ level. The results of radiative transfer computations 
(Hydrolight code) show that a maximum difference of 5% 
exists between K^ 9 and K 9 f' for a TChl a concentration 
of 0.3 mgm -3 and a solar zenith angle of 45°, which trans- 
lates into a maximum error in L u ( 0 _ ) of 1.25% when using 
K^’ 9 t° perform the extrapolation, and when Kp is about 
0.05m -1 (i.e., a TChl a concentration of 0.3 mgm -3 ). This 


f The distance between the two depths of measurements is 
exactly known; the depth in question here is the depth Z 
from which the extrapolation to just below the surface is 
performed through L u ( 0 _ ) = L u (Z)e KljZ . 


error could be significantly reduced by correcting the val- 
ues of K°£ 9 following what can be learned from the results 
of radiative transfer computations. 

At the end, assuming an average error of 3% on the 
estimation of L u ( 0 _ ), because of uncertainties in the ex- 
trapolation to the 0“ level, seems realistic. 

10.2.1.2 Bidirectional and Transmission Effects 

The in-water measurements provide the upwelling ra- 
diance at nadir, which has to be transformed into the up- 
welling radiance for the direction below water (O') corre- 
sponding to an above- water zenith angle ( 0 ), itself corre- 
sponding to a given viewing angle from the satellite, 0 V . 
The transformation is simply performed by ratioing the 
Q-factor at nadir to the Q-factor for the direction O'. The 
uncertainty is here only in the relative values of Q at these 
two directions (in a match-up configuration, the difference 
between both will be minimized, i.e., 6 < 20°). In the 
blue domain, where the geometry of the light held does 
not depend very much on the particle phase function, the 
ratio of the Q-factors is probably correct to within a small 
amount, so a 2% uncertainty is assumed here. 

The last step to get the water-leaving radiance is to 
multiply L u (0') by the expression (1 — p(6'))n~ 2 , which 
accounts for the transmission across the interface, where 
p(0') is the Fresnel reflection coefficient for the water-air 
interface and n w is the refractive index of water. This step 
does not introduce any error as far as O' < 20° and the sea 
surface is approximately hat (wind speed less than 15kts). 

10.2.1.3 Atmospheric Transmittance 

The water-leaving radiance, Lw, has to be multiplied 
by the atmospheric diffuse transmittance, t,p before it can 
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be added to the atmospheric path radiance. If the ex- 
pression in use for computing t, ( j is considered valid (Wang 
1999), the largest source of uncertainty is the aerosol op- 
tical thickness, which has a small impact on t,i. A factor 
of two change in r a (from 0.1 to 0.2), for example, trans- 
lates as a 2% error on td at 443 nrn, which results in a 2% 
uncertainty on the product tdLw 

In the green domain (e.g., at approximately 550 nrn), 
another source of uncertainty is the ozone content of the 
atmosphere. When A = 550 nm, the Rayleigh optical thick- 
ness is about 0.09 and the ozone optical thickness is ap- 
proximately 0.03 for an ozone content of 350 DU. The un- 
certainty in td{ 550) corresponding to an uncertainty of 
50 DU would be of the order of 0.5% for viewing angles 
less than 45° (with, for instance, td = e~ ( '°' 5Tr+Toz ^ cos9v ). 


AOT is not zero. A 0.2% uncertainty is assumed here when 
the atmospheric pressure is known (local measurements 
from the Meteo buoy near the BOUSSOLE site). 

10.2.3.2 Sea Surface State 

Quantification of the uncertainty due to an incorrect 
representation of the sea surface state is not really possi- 
ble, and it would be probably useless because of the large 
uncertainties that remain in the parameterization of the 
surface wave slope probability distribution as a function of 
the wind speed. A tentative value of 0.5% is assumed here 
for the uncertainty introduced by these surface effects in 
the computation of L t at the TOA level, with the under- 
lying assumption that experiments are only performed for 
days of calm weather. 


10.2.2 Above- Water Measurements 


10.2.3.3 Aerosol Optical Thickness 


The uncertainty to which L\ y is derived from the radi- 
ance measured above the sea surface (which includes L\y 
and surface reflection effects), can be considered as be- 
ing on the order of about 10% (Hooker et al. 2002b and 
Hooker and Morel 2003), even when performed in well- 
controlled experimental conditions and in excellent envi- 
ronmental conditions. A careful implementation of the re- 
quired metrology and proper removal of glint effects during 
data processing can reduce the uncertainties to the same 
order as the calibration uncertainty (Hooker et al. 2004) 
for above- or in-water radiometers. 

Because above-water measurements are collected well 
above the sea surface, the uncertainty attached to the ex- 
trapolation of in- water measurements from depth to just 
beneath the surface does not enter into account. The 
same uncertainty is attached to the determination of bi- 
directional effects ( Q factors), as well as to the calculation 
of the atmospheric diffuse transmittance. 


10.2.3 Atmospheric Path Radiance 

The uncertainties assoicated with computing the at- 
mospheric path radiance are examined based on the un- 
certainties in the parameters used in the computations. 


10.2.3.1 Atmospheric Pressure 

Atmospheric pressure, P, is known to within 0.5% (i.e., 
5hPa). The changes in the path reflectance due to changes 
in P can be expressed as (Antoine and Morel 1998): 


ATpath 

A P 


1 + 


A P 

j> 'lr. 


(17) 


where ly is the ratio of the Rayleigh optical thickness to the 
total (aerosol plus Rayleigh) optical thickness. In the worst 
case, where r) r = 1 (no aerosols), the uncertainty in P is 
transferred onto T pa th) but is actually a little less as shown 
by Gordon et al. (1988). It will also be less as soon as the 


The AOT should be known to less than 0.01 in absolute 
units (Fargion and Mueller 2000). The relative changes in 
Lpath due to changes in AOT in the blue and green domains 
(i.e., around 440 and 550nm, respectively), AL pa th/Ar a , 
are of the order of 0.1, e.g., as shown in Fig. 5 in Antoine 
and Morel (1998). An uncertainty of 0.01 in r a should, 
therefore, translate as an uncertainty of 0.001 in L pa th) 
which corresponds to about 1% of L pa th(4:43) and 2% of 
L path (550) when r a = 0.1. The uncertainties are slightly 
less if r a is larger. These two values (1% and 2%) are 
adopted here, while noting that the corresponding numbers 
become 2-3% at 865 nm. 

10.2.3.4 Aerosol Type 

Retrieving the aerosol type, actually the aerosol volume 
scattering function (VSF) and single scattering albedo, is 
much more difficult than retrieving the AOT. The reasons 
for the difficulty are, a) the retrieval is based on inversion of 
sky radiance measurements, which are delicate to perform 
with the desired calibration constraints, and b) the inver- 
sion itself necessarily uses assumptions and is also subject 
to uncertainties. 

In order to minimize the impact of uncertainties on 
the aerosol type, situations where the satellite data and 
external information (wind direction, independent obser- 
vations, etc.) are both available are selected, in such a 
way that there is a high likelehood of the presence of mar- 
itime aerosols. Note that information about the relative 
humidity (RH) at the sea level will be available from the 
measurements of the Meteo buoy located in the vicinity of 
the BOUSSOLE mooring. 

It is, therefore, considered here that the uncertainty is 
entirely due to the uncertainty on RH (the value itself or 
the vertical profile), for maritime aerosols conforming to 
the description of Shettle and Fenn (1979). From the re- 
sults of radiative transfer simulations performed with these 
models, it appears that an uncertainty of 15% in the rel- 
ative humidity (RH equal to 70% or 99% instead of 80%) 
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leads to a 2% uncertainty in the computation of L pat h(443) 
when r a < 0.1 and 9 S < 40°. This uncertainty is closer to 
5% in the NIR (A = 865 nrn). 

Maintaining this level of uncertainty means that the 
vicarious calibration experiments have to be performed for 
T a values as low as possible and for high solar elevations 
(the minimum solar zenith angle at the BOUSSOLE site 
is about 21°). 

10.2.3.5 Ozone Content 

The importance of ozone content is investigated here 
in terms of its effect on the calculation of T pa th (the other 
effect has been already considered in the calculation of 
the atmospheric diffuse transmittance). Ozone does not 
strongly affect the blue bands, which are the most critical 
for ocean color applications. A null impact on L pa th is 
considered here. 

For the green band, the problem is different because 
the ozone absorption is there at its maximum, with an 
optical thickness of about 0.03 at 550 nm when the total 
ozone amount is equal to 350 DU. The impact on L pat h 
remains, however, extremely weak, i.e. , less than a 5% 
difference between L pat h for an atmosphere with a stan- 
dard ozone amount (350 DU) and an ozone-free atmos- 
phere. The typical uncertainties in the ozone amount (i.e., 
±20 DU) should, therefore, be without significant impact 
on the computation of L pat h (less than 0.5%). 

The calibration uncertainty has been already consid- 
ered when examining the uncertainty due to incorrect de- 
termination of the AOT and aerosol type. In principle, a 
valid radiative transfer code, when fed with accurate input, 
provides an exact answer in terms of radiance. It is, there- 
fore, assumed that the modeling itself does not introduce 
any additional uncertainty. 

10.2.4 Marine Contribution to Total Radiance 

In clear oligotrophic waters, L\y(\) reaches a maximum 
value in the blue part of the spectral domain, and corre- 
sponds to about 15% of the TOA total radiance at 440 nm 
when the AOT is 0.05. Conversely, Lw(\) experiences 
minimum values in the green domain, and corresponds, 
for example, to about 5% of the TOA total radiance at 
550 nm when the AOT is 0.3. Consequently, the various 
uncertainties in the estimation of L\y( A) as given above, in 
terms of percentages, correspond to uncertainties in L t ( A) 
ranging from 0.15% to 0.05%. 

The various terms are summarized in the Tables 12 and 
13, and are based on three different configurations concern- 
ing the estimation of the marine and atmospheric contribu- 
tions to the TOA total radiance: in-water measurements of 
L u ( A) and measurements of the aerosol properties, above- 
water measurements of Lw(X) and measurements of the 
aerosol properties, and, finally (and tentatively), estimated 
values for the water-leaving radiance and aerosol proper- 
ties (the so-called Rayleigh scattering above clear water 
technique). 


10.2.5 Conclusions 

In summary, the values in Tables 12 and 13 show that 
a 2-3% uncertainty in the simulation of the TOA total 
radiance can be reached in extremely favorable conditions. 
Considering the possible type and range of uncertainties 
that are possible in the entire exercise, a more realistic 
estimate of the accuracy to which the vicarious calibration 
can be performed would be approximately 5% in the visible 
part of the spectrum. 

10.3 A Practical Example 

The inversion algorithm suggested in Sect. 7.7 not only 
permits the characterization of the microphysical prop- 
erties of the aerosols providing the refractive index, but 
also supplies knowledge about the aerosol phase function, 
which is a key parameter in the estimation of the atmos- 
pheric radiance and transmittance. The proposed algo- 
rithm can thus be used to lead a vicarious calibration of 
ocean color sensors. 

The focus here is on the vicarious calibration of the 
MERIS NIR band at 865 nm. The sun photometer (Sects. 
3.3, 4.3, and 5.3) is used to provide the necessary informa- 
tion about the aerosol properties. The procedure for the 
vicarious calibration is outlined as follows: 

1. The aerosol refractive index is derived with the in- 
version algorithm. 

2. The optical depth and Angstrom exponent are mea- 
sured. 

3. The TOA radiance is computed with the OSOA 
radiative transfer model for similar geometries of 
MERIS scenes on the basis of the aerosol optical 
properties, which is compared with the TOA signal 
recorded by MERIS. 

Some restrictions and assumptions should be made, how- 
ever, to apply the vicarious calibration technique: 

4. The ocean is assumed to be black at 865 nm which 
is realistic in the Mediterranean Sea (open-ocean 
waters), so no marine signal is included in the com- 
putations. 

5. The MERIS images are selected for geometrical con- 
ditions based on avoiding glint in the vicinity of the 
ground photometer. 

6. Because of the proximity of the land in each im- 
age, the aerosol properties may exhibit a strong 
spatial variability at the land-ocean interface, so 
the aerosol model derived from the coastal sun pho- 
tometer will be applicable to MERIS marine pixels 
only if the aerosols in the study area are spatially 
homogeneous (this constraint substantially reduced 
the number of available images for the calibration 
exercise) . 
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Table 14. Computational aspects used in the OSOA radiative transfer model for the calculation of the top of 
atmosphere radiance in the MERIS geometry: r a at 675 nm (measured); a (measured); 9 S (in degrees); viewing 
zenith angle (in degrees); relative azimuth, A cj), between sun and satellite (in degrees); scattering angle, Od , 
(in degrees); aerosol model (complex refractive index m r — jrrii) derived by the neural network. The time 
corresponding to MERIS overpass (t M ) and to the chosen sequence of principal plane measurements ( t p ) is 
indicated in GMT (the year is 2002). 


Date 


t P 

To, 

a 

0s 

o v 

A cj) 

o d 

m r 

- jrm 

11 July 

1000 

0839 

0.134 

1.34 

29.5 

7.7 

25.3 

157.2 

1.46 

- j'0.009 

7 September 

0935 

1616 

0.097 

1.59 

45.1 

39.0 

38.6 

153.8 

1.41 

- jO.OOO 

26 September 

0938 

1027 

0.041 

1.29 

50.5 

35.7 

45.4 

146.4 

1.55 

- J0.004 

1 October 

1021 

1125 

0.074 

1.61 

48.7 

26.8 

123.3 

113.8 

1.52 

- jO.Oll 

12 October 

1021 

1122 

0.037 

1.02 

56.0 

39.1 

49.4 

140.7 

1.57 

- j'0.008 


7. For each selected image, a subscene of 10x10 pixels 
is extracted far enough from the coastline (typically 
at a distance greater than 6 km) to avoid the adja- 
cency effects on the radiance measured by the satel- 
lite sensor caused by the radiometrically brighter 
land (Santer and Schmechtig 2000). 

8. After correction for the eccentricity of the Earth 
orbit, the TOA radiance is averaged over this sub- 
scene. 

The inversion algorithm is applied on ground-based 
measurements collected at the time of the satellite over- 
pass (around 1000 GMT for the Villefranche site) so the 
simulated radiance is consistent with MERIS data. Never- 
theless, occasionally measurements are not available at the 
time of overpass; either they are missing or of bad qual- 
ity. In such situations, the aerosol model derived from the 
inversion of the closest (in time) measurements of satellite 
overpass is used to reconstruct the TOA signal. This is 
possible as long as stable conditions are observed during 
the day. The daily variation of the Angstrom exponent is 
used to check this assumption. 

10.4 Preliminary Results 

Table 14 lists the computational aspects used in the 
OSOA model (geometry, aerosol refractive index, and op- 
tical depth) to predict the TOA radiance corresponding 
to MERIS observations. Only five matchups were avail- 
able at the time of this preliminary analysis, depending on 
the quality of ground measurements and MERIS scenes. 
It should be noted that on 7 September 2002, the ground 
measurements made on the morning showed questionable 
data. Because r a (675) and a were stable all along this 
day — r a (675) varied from 0.097 to 0.119, and a varied from 
1.45 to 1.69 — the aerosol model derived from the measure- 
ments collected at the end of the day (at 1616 GMT) was 
used, because it includes the maximum value of the scat- 
tering angle (150°). 

Table 15 presents a comparison between MERIS TOA 
radiances and OSOA calculations. The agreement is very 


good and remains below 5%. Gordon (1998) showed that 
for a given absolute calibration uncertainty in the NIR 
domain, the calibration error progressively decreases with 
decreasing wavelength. More recently, Wang and Gor- 
don (2002) demonstrated that a calibration uncertainty of 
±15% at 865 nm leads to an acceptable calibration uncer- 
tainty smaller than 3% at 412 nm. The results presented 
here are, therefore, within the accuracy requirements for 
the NIR band. Those results are still preliminary, how- 
ever, because the number of available match-up data is 
not significant enough to draw rigorous conclusions about 
the MERIS calibration. 

Table 15. The comparisons between the TOA 
radiances computed from OSOA radiative transfer 
calculations (±o) or from MERIS (Lm) at 865 nm 
(expressed in units of mWm“ 2 nm" 1 sr" 1 ). The 
relative difference, AL/L, corresponds to the ratio 
( Lm — Lq)/Lm- 


Date 

Lm 

Lo 

A L/L [%] 

11 Jul. 

4.4248 

4.2194 

4.64 

7 Sep. 

4.2832 

4.4165 

-3.11 

26 Sep. 

2.7727 

2.8420 

-2.50 

1 Oct. 

2.7479 

2.7362 

0.42 

12 Oct. 

2.7550 

2.7337 

0.77 


As a future activity, the vicarious calibration method 
will be extended to the visible bands. In the visible part of 
the spectrum, the marine signal is not equal to zero, so the 
additional difficulty will be to get simultaneous in-water ra- 
diometric measurements in order to reconstruct the total 
TOA radiances as the sum of the atmospheric signal (cal- 
culated using the same technique as the NIR band) and 
the marine signal multiplied by a diffuse transmittance. 
To achieve this task, the continuous record collected by 
the BOUSSOLE buoy will provide the needed information 
on the signal exiting the water and will increase the num- 
ber of matchups needed to perform a relevant vicarious 
calibration of MERIS. 
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Glossary 


AC-9 

AC-9+ 


AERONET 

AOP 

AOPEX 

AOT 

ARGOS 


ASCII 


Absorption and Attenuation Meter 
Absorption and Attenuation Meter (with data 
handling) 

Aerosol Robotic Network 

Apparent Optical Property 

Advanced Optical Properties Experiment 

Aerosol Optical Thickness 

Not an acronym, but the name given to the 

data collection and location system on the 

NOAA operational satellites. 

American Standard Code for Information In- 
terchange 


BOUSSOLE Bouee pour Tacquisition de Series Optiques a 
Long Terme (buoy for the acquisition of a long- 
term optical series). 


C-star 

CDOM 

CE-317 

CE-318 

CERT 

CIMEL 

CNES 

CNRS 

Coasts 

ComPACT 

CTD 

CZCS 


WETLabs transmissometer 

Colored Dissolved Organic Matter 

CIMEL Electronique (sun photometer model) 

317 

CIMEL Electronique (sun photometer model) 

318 

Calibration Evaluation and Radiometric Test- 
ing 

Not an acronym, but the name of a French com- 
pany manufacturing sun photometers. 

Centre National des Etudes Spatiales (French 
space agency). 

Centre National de la Recherche Scientifique 
(French national center for scientific research). 
Costal Atmosphere and Sea Time Series 
Compact Portable Advanced Characterization 
Tank 

Conductivity, Temperature, and Depth 
Coastal Zone Color Scanner 


DACNet 
DAS 
DATA- 100 
DCM 

DHI 

DOC 

DP 

DUT 

DYFAMED 


ECO-BB3 

ESA 

EZ-III 

FEL 

GAC 

GF/F 

GMT 

GPS 

GUI 

HOBI Labs 

HPLC 

Hydroscat-II 


Data Aquisition and Control Network 
Data Acquisition Sequence 
(Satlantic) Data (acquisition series) 100 (unit) 
Deep Chlorophyll Maximum 
Not an acronum, but the name of the DHI Wa- 
ter and Environment Institute in Denmark. 
Dissolved Organic Carbon 
Diagnostic Pigment 
Device Under Test 

Dynamique des Flux Atmospheriques en Medi- 
terranee (dynamics of atmospheric fluxes in the 
Mediterranean) . 

Three-channel backscattering (meter) 
European Space Agency 
Compass and two-axis tilt sensor 

Not an acronym, but a lamp designator. 

Global Area Coverage 

Not an acronym, but a specific type of glass 
fiber filter manufactured by Whatman. 
Greenwich Mean Time 
Global Positioning System 
Graphical User Interface 

Hydro-Optics, Biology, and Instrumentation 
Laboratories 

High Performance Liquid Chromatography 
HOBI Labs Hydroscat-II backscattering meter 


ID 

IFREMER 


INSU 


IOCCG 

IOP 

IR 

JGOFS 

L19 

LCD 

LED 

LOV 


Identification 

Institut Frangais de Recherce pour I’Exploita- 
tion de la Mer (French research institute for 
exploitation of the sea) 

Institut National des Sciences de TUnivers 
(French national institute of sciences of the uni- 
verse). 

International Ocean Color Coordinating Group 

Inherent Optical Property 

Infrared 

Joint Global Ocean Flux Study 

(Perkin Elmer) Lambda 19 (spectrometer) 
Liquid Crystal Display 
Light-Emitting Diode 

Laboratoire d’Oceanographie de Villefranche 
(Oceanographic Laboratory of Villefranche) 
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MARINTEK Marine Technology Research Institute 

MEGS MERIS Ground Segment Simulator 

MERIS Medium Resolution Imaging Spectrometer 
MINItracka Not an acronym, but a fluorometer manufac- 
tured by Chelsea Company. 

MLAC Merged Local Area Coverage 

MOBY Marine Optical Buoy 

MODIS Moderate Resolution Imaging Spectroradiome- 
ter 

MODIS - A MODIS instrument on the Aqua satellite. 

MVDS Multichannel Visible Detector System 

NASA National Aeronautics and Space Administra- 
tion 

NIR Near-Infrared 

NIST National Institute of Standards and Technol- 
ogy 

NOAA National Oceanic and Atmospheric Adminis- 
tration 

OCL200 Satlantic Ocean Color Irradiance-200 series 
(sensor) 

OCTIOOO Satlantic Ocean Color Irradiance-1000 series 
(sensor) 

OCR-200 Satlantic Ocean Color Radiance-200 series (sen- 
sor) 

OSOA Ordres Successifs Ocean Atmosphere (ocean at- 
mosphere successive orders) 

PAR Photosynthetically Available Radiation 
PARASOL Polarization and Anisotropy of Reflectances for 
Atmospheric Sciences coupled with Observa- 
tions from a Lidar 
PC Personal Computer 

POLDER Polarization and Directionality of the Earth Re- 
flectance 

PROSOPE Productivite des Systemes Oceaniques Pelag- 
iques (productivity of pelagic oceanic systems) 

R/V Research Vessel 
RAM Random Access Memory 
RH Relative Humidity 
rms Root Mean Square 

SBE Seabird Electronics 

SCOR Scientific Committee for Oceanographic Re- 
search 

SeaHARRE SeaWiFS HPLC Analysis Round-Robin Exper- 
iment 

SeaWiFS Sea-viewing Wide Field-of-view Sensor 

SIMBAD Satellite Validation for Marine Biology and 
Aerosol Determination 
SIMBADA Second-generation SIMBAD 

SIMBIOS Sensor InterComparison and Merger for Bio- 
geochemical Studies 

SIRREX SeaWiFS Intercalibration Round-Robin Exper- 
iment 

SMSR SeaWiFS Multichannel Surface Reference 

SPMR SeaWiFS Profiling Multichannel Radiometer 
SQM SeaWiFS Quality Monitor 

SQM-II Second generation SQM 
SST Sea Surface Temperature 

STAM Satlantic Telemetry Acquisition Manager 

TOA Top of Atmosphere 
TSM Total Suspended Matter 


UPMC Universite Pierre et Marie Curie (translated as 
University of Pierre and Marie Curie) 

UV Ultraviolet 

VSF Volume Scattering Function 

WCF WETLabs CDOM Fluorometer 
WETLabs Western Environmental Technology Laborato- 
ries 

Symbols 

a( A) The spectral absorption coefficient. 
ad( A) The spectral absorption coefficient for detritus. 
a p ( A) The spectral absorption coefficient for all particles. 
a^(A) The spectral absorption coefficient for phytoplank- 
ton. 

6(A) The spectral scattering coefficient. 

6b (A) The spectral backscattering coefficient. 

c(A) The spectral attenuation coefficient. 

C a The chlorophyll a concentration. 

Ed{ A) The spectral downward irradiance. 

FA, (A) The spectral downward irradiance above the sea 
surface. 

E u ( A) The spectral upward irradiance. 

/ The /-factor ( R/{bb/a )). 

Kd{ A) The spectral diffuse attenuation coefficient for 
downward irradiance. 

Kl{ A) The spectral diffuse attenuation coefficient for up- 
welling radiance. 

K u ( A) The spectral diffuse attenuation coefficient for up- 
ward irradiance. 

L g The radiance due to sun glint at the TOA. 

Fpath(A) The spectral path radiance at the TOA. 

Lt The total radiance at the TOA. 

L u ( A) The spectral upwelling radiance. 

Lw{ A) The spectral water- leaving radiance. 

[Lvv(A)] n The spectral normalized water- leaving radiance. 
L wc The radiance due to white caps at the TOA. 

m The slope of a linear regression fit. 
mi The imaginary part of the complex refractive index 
of aerosols. 

m r The real part of the complex refractive index of 
aerosols. 

n w The refractive index of sea water. 

P The atmospheric pressure. 

Q( A) The spectral Q-factor ( E u /L u ). 

Q n The Q-factor at nadir (O' = 0°). 

f The average reflection coefficient for upward irradi- 
ance from below the sea surface. 
r 2 The regression coefficient of determination. 

R The irradiance reflectance. 
ik(0') A geometrical factor accounting for all refraction 
and reflection effects at the air-sea interface. 

3?o St(0' = 0). 

rd The radius of the SPMR radiometer. 
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T The total transmittance along a pixel-to-satellite 
sensor path. 

td The diffuse transmittance along a pixel-to-satellite 
sensor path. 

[TChl a] The total chlorophyll a concentration. 

y 0 The y-intercept of a linear regression fit. 

2 Depth. 

a The aerosol Angstrom exponent. 

At/) The relative azimuth difference angle. 

Tf r The ratio of the Rayleigh optical thickness to the 
total (i.e., aerosol plus Rayleigh) optical thickness. 

O' The in-water refracted viewing angle sin -1 0 V / 1.34. 

Od The scattering angle. 

Os The solar zenith angle. 

0 V The satellite viewing zenith angle. 

A Wavelength. 

p The average reflection coefficient for downward ir- 
radiance across the air-sea interface. 
p(0) The Fresnel reflection coefficient for incident angle 
0 . 

T a Aerosol optical thickness. 
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